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1, GENERALIDADES 


1.1, INTRODUÇÃO 


O tipo de problema de interesse da dinâmica estrutural & 
o que incorpora modificações nas quantidades de movimento dos 
sistemas elasticos. Assim, a ação de um motor sobre sua base, a 
do vento ou das ondas do mar, a de terremotos ou explosões,a de 
impacto ou de cargas moveis, sobre sistemas estruturais, são 
exemplos. desses casos, em que o equilibrio dos sistemas so é ve 


. E Se (v'auce Her) 
rificado com a consideração de forças de massa. 


A variação com o tempo das forças que atuam sobre um sis 
tema deformavel, faz com que, não sendo o desenvolvimento das 
forças elasticas suficientemente rapido para manter o equili- 
brio, o sistema modifique a sua situação cinemática para buscar 
o equilibrio com a ajuda das "forças de massa", ou “forças de 
incrcia". Em linguagem mais precisa, o sistema tem sua quantida 
de de movimento alterada, e a 23 lei de Newton assegura a satis 


fação dus condições de equilibrio. 


Basicamente, a grande modificação e a necessidade de in- 


corporação da variavel tempo, t, às equações de equilibrio. 


0 movimento do sistema transforma-se numa oscilação, pe- 
la sucessiva troca de energia potencial em cinetica, e vice-ver 
sa, e & dito que a estrutura “vibra”. Alem do carater repetiti- 
vo, ou alternativo, da resposta, as amplitudes dos deslocamen- 
tos chegam a ultrapassar, em diversas vezes, os valores corres- 
pondentes à aplicação estática da ação, o que empresta ao estu- 


do de dinâmica estrutural uma importancia indiscutível. 


Toma-se, por exemplo, a viga representada na fig. 1.1;na 


parte (a) dessa ilustração encontra-se o gráfico do deslocamen 


to vertical do ponto c, para uma força A agindo estaticamente; 


na parte (b) da figura, mostra-se o diagrama do mesmo desloca- 


mento para uma força subita de mesma intensidade Fa 


F(t) + Nai 
Ve 
A c | 
Pecas de 
42 Teen icoríta 
F(t) 
Fo 
y t 
vit)/ve 


Fig. 1.1. Respostas de cargas: (a) estática; (b) dinâmica 


Verifica-se que a simples transformação da carga estati- 
ca numa ação subita faz o deslocamento do ponto c ser repeti- 
do e com uma amplitude duas vezes maior que no caso estatico. 
Compreende-se que o exemplo considerado e o de um sistema idea] 
onde não ocorre dissipação de energia - nas estruturas reais 
ha elementos, ou reações internas, que promovem um consumo de 
energia que provoca uma redução, com o tempo, quer da amplitu- 


de, quer do carater oscilante do movimento. 


Um outro ponto peculiar são as acelerações a que o siste 


ma fica submetido, durante a vibração, e que podem limitar a sua 


ocupação, comoda ou estavel, por pessoas, itens de equipamento, 


etc. 


1.2. CLASSIFICAÇÃO 


As ações dinamicas que ocorrem em um sistema estrutural 
podem ser avaliadas de duas maneiras: detenministica, na qual as 
caracteristicas do sistema e da excitação são à priori estabele 
cidas e não-deterministica, randômica ou estocastica quando pe- 
lo menos um desses elementos é estabelecido a partir de um con- 


junto amostral com certa definição probabilística. 


A vibração € um termo que descreve a oscilação num siste 
ma mecânico, e na pratica não possui muitas vezes um padrao 
regular, podendo ser uma combinação de vários harmônicos de res 
posta simples. Se ela repete-se a certos intervalos de tempo & 


dita periódica, do contrário E não-períúdica, ou complexa. 


Muitas vezes, a vibração é um movimento periódico, no qual 
todas as suas características se repetem num determinado inter 
valo de tempo denominado periodo. A fig. 1.2 mostra dois movi - 
mentos periódicos, qualquer (a) e harmônico (b); neste último ca- 
so, pode ser escrito o seguinte: v=v, cos w t, onde ug=2n ta e 

ad 1 €3 lar, f.s L e a frequencia nat 1 
a frequencia natural circu + Ta” T q ura e 


To, & o periodo natural. 


As ações atuantes em um sistema material podem ser tam 
bêm, de curta dunação, na qual as forças dinamicas externas ao 
sistema agem num pequeno intervalo de tempo, isto e, numa fra- 


ção do período natural, caso contrário são de Tonga duração. 


As vibrações podem também ser classificadas em Livres, 


quando provocadas exclusivamente pela energia potencial e cine- 


tica existentes no sistema, ou tourçadas, que incorporam ainda as 


ações externas sobre o sistema, variaveis com o tempo;e ainda,a 


vibração é amortecida quando o sistema dissipa energia durante 


a oscilação, sendo nao-amonrtecida no caso oposto. E importante 


salientar que todo sistema nao-amortecido é uma idealização pa- 


ra auxiliar no entendimento de situações limites do caso amorte 


cido. 


f(4) (1) 


— 


(o) 


Fig. 1.2. Movimentos periódicos: (a) qualquer; (b) harmônico 


2, SISTEMAS COM UM GRAU DE LIBERDADE (SIGL) 
2.1. CONCEITUAÇÃO 


2.1.1. DEFINIÇÃO 


Pode-se definir grau de liberdade dinâmico de um sistema 
Ad MINAMISO 


E a Eis 


da seguinte maneira: "numero minimo de coordenadas para se saber 


— 


como se encontra, em qualquer instante, a configuração do siste- 


ma". Então para SIGL deve-se ter uma única coordenada. 
| | 
a I 
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Fig. 2.1. Exemplos de sistemas dinamicos: (a) sistema de massa 
simples; (b) pêndulo simples; (c) viga de rigidez infinita com 
apoio elástico a direita; (d) disco rigido sujeito à torção;(e) 
viga de massa desprezível com duas massas puntiformes, m, e m,; 
(f) viga de massa desprezível com uma massa não-puntiforme; (g) 
viga com massa distribuida ao longo da coordenada x 


Na fig. 2.1, & mostrada uma série de exemplos de sistemas 
dinâmicos. Na parte (a), o sistema estã sujeito somente a deslo 
camentos verticais, pois a massa é obrigada a se deslocar na di- 
reção y, possuindo 16G.L. Caso não houvesse restrições das pare- 
des, haveria 6G.L. e ainda no caso da massa da mola não ser des- 
prezível, o numero de G.L. seria infinito - & sô pensar que cada 
ponto possuiria 6G.L. Outros exemplos de SIGL podem ser vistos 
em (b), (c), (d). Nas partes (e) e (f) tem-se S2GL e em (g) um 
sistema com um numero infinito de G.L.- neste caso se & conhecida 


a configuração na qual a viga ira vibrar, o sistema tem 1G.L.. 
2.2. EQUAÇÃO DO MOVIMENTO 


2.2.1, O OSCILADOR SIMPLES 


0 oscilador simples e um modelo de estudo para SIGL e & 
formado de massa, mula e amortecedor. A fig. 2.2a mostra este sis 
tema onde k & a constante elastica, ou rigidez da mola-força ne- 
cessiria para um deslocamento unitário; (c/& o coeficiente de amor- 
tecimento viscoso - força correspondente a uma velocidade unita- 


ria. As hipoteses para este modelo são as seguintes: 


1) a mola tem massa desprezível; “ 
2) as resistênciasoferecidaspela mola e pelo amortecedor são 
proporcionais ao deslocamento e à velocidade, respectiva- 
PF; - 0 Lu > diadR 
mente; ” cr lmwioudade) 


3) não hã perdas de energia devido a atrito que atua externa 


mente ao sistema. 


2.2.2. MÉTODOS DE ANÁLISE 


2.2.2.1. EQUILÍBRIO DIRETO 


Na parte (b) da fig. 2.2, sao mostradas as forças de iner 


cia, da mola, do amortecimento e a externa, respectivamente ro 
E» Fa F(t). A equação de movimento para o sistema pode ser ob 
tida a partir da 2a lei de Newton ou pelo princípio de D'Alem- 


F 


bert, como se segue. 


mu(t) + cult) + kv(t) = F(t) EB,1) 


A equação obtida e diferencial, ordinaria, de segunda or 


dem, linear, não-homogênea e com coeficientes constantes. 


A fig. 2.3 mostra o mesmo oscilador, com as forças atuan 
do na direção vertical, parte (b), com a inclusão do seu peso,pP, 
representando todo um conjunto estatico de cargas. Assim sendo, 


a equação diferencial obtida e a seguinte: 


mu(t) + cult) + k(v(t)+v.c) = F(t) + P (2.2) 


est 


Na expressão acima, vAcy= P/k e o deslocamento estatico 


es 


do sistema. Substituindo-se esta definição na equação (2.2) ob- 


tem-se a mesma equação (2.1). 


klvevest) MV ey 


* Nest POSIÇÃO DE EQUILIBRIO 
ESTATICO 


+, 
, 
Les Lo-—s 


P+ F(t) dis 


(a) (b) 


Fig. 2.3. Oscilador simples com a consideração do peso 


Isto significa que, ao estudar-se a equação dinâmica, po 
de-se considerar as forças relativas a partir da posição de equi 
Níbrio estático. As forças estaticas, como visto acima, podem 
ser desconsideradas, pois cancelam-se mutuamente. Deve-se ter em 
mente, todavia, que para fins de dimensionamento, a força da mo 


la deve ser considerada, utilizando a expressão seguinte: 


F=ky(t)+ Po. (2.3) 


2,.2,2,.2, PRINCÍPIO DOS TRABALHOS VIRTUAIS 


No modelo apresentado na fig. 2.2a, a massa agora sofre 
um deslocamento virtual, nao nulo, e dessa forma as forças rea- 


lizam trabalho. 


O trabalho total, virtual, realizado pelo sistema deve 
ser nulo, ou seja: 


ana 6v - F, 6v - Fe 6v + F 6v =0 (2.4) 


A 


(-Fy-Pa-Fp4FDav = O 


O sinal negativo significa que a força tem sentido contra 


rio ao deslocamento virtual. 


O resultado encontrado & a mesma equação (2.1). 


2,2,2.3, PRINCÍPIO DE HAMILTON 


O princípio de Hamilton pode ser expresso pela seguinte 
formulação variacional: 
Ei t, 
E e(T-V)dt + é Vc dt = 0, (845) 


ts t, 
considerando-se o mesmo oscilador da fig. 2.2, tem-se: 
T= | mv (2,541) 
7) > 
energia cinetica do sistema; 
tela , kv? (2.5.2) 


energia de deformação da mola, V, que e igual à energia potencial do 


sistema, U; 
uno = Flt)6v - cv 6v (2.6.8) 


trabalho virtual realizado pelas forças não conservativas. 


Substituindo-se as expressões (2.5.1), (2.5.2), (2.5.3) na 
expressão (2.5), obtem-se: 


t, 
| [mv. 6v-cv. 6v-kv. Gv+F(t)sv]dt = 0 (2.6) 
t 


] 


Integrando-se, por partes, o primeiro termo da eg. (2.6), 
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notando-se que 6v=d(6v)/dt, vem: 


A” o Es a is 
J mv 6v dt = mv &v do | my 6v dt (2.7) 
E t t 


Desde que seja adotado que, no princípio de Hamilton,sv 
os E 2 

desaparece nos limites de integração, o termo mvov|, se anula 
1 


e assim: 


j my 6v dt = -| mv 6v dt (2.8) 
t T 


Voltando, com a exp. (2.8), à eq. (2.6), tem-se: 


t, 
| (-my-cv-kv+F(t))óv dt = 0 (2.9) 
t 


1 


Como a variação 6v & arbitrária, pode ser eliminada da in 


tegral e fica claro que a equação de movimento & a mesma da ex 


pressão (2.1). 


2.2.3. OUTROS EXEMPLOS PARA SIGL 


Exemplo 2.1. É mostrada, na fig. 2.4, uma viga de massa despre- 
zivel, rígida, de vão L, sujeita a uma força F(t), aplicada a 
uma distância D do apoio rotulado O'. A mesma possui um apoio 
xTEKich é un amortecedor distanciados, respectivamente, A e B 
de O' e uma massa concentrada na extremidade direita, como po- 
de-se ver em (a). A parte (b) da ilustração mostra todas as 
forças agindo sobre o sistema na posição de equilibrio estati- 
co. Im (c), os deslocamentos dinâmicos relativos à posição de 


equilíbrio estático sao apresentados e em (d) aparecem as for- 


1] 


ças dinamicas atuantes na viga. 


VIGA RÍGIDA S/ MASSA Ftt) 
o m 
= k he : P + P. mg 
; - A 


, o EQUILÍBRIO ESTÁTICO 
B 


(c) (d) 


Fig. 2.4. Exemplo 2.1 


A equação de movimento pode ser obtida da forma seguinte, 


* utilizando-se do princípio de D'Alembert, isto &, fazendo-se o 


somatório dos momentos em relação a 0' nulo: 


: . B2 A? 
; E > mLy cv + kv5>—-= F(t)D (2.10) 

EM O + Te L 

2 2 

co mir eb) v + PD) v = 7 FCt) (2.11) 
my + Ceq Vv + keq V = Feq (2.12) 
Ba” ha” .D - : 
onde Ceg * c(T) ; keq = k(p) e Feq SE F(t) são os valores equi 


valentes de amortecimento, rigidez e força externa; respectivamente. 


A partir da eq. (2.10), utilizando-se a transformação 


v=L6, segue-se: 
” as A,* D 
mó + c(r) 8+ k(D) 0 = Tr F(t) (2. 14) 


Exemplo 2.2. A mesma viga do exemplo 2.1, é colocada agora na 


posiçao vertical, fig. 2.5, considerando-se uma força P atuando 


verticalmente. 


VIGA RÍGIDA 
SEM MASSA 


o | EQUILISMO| ESTÁNCO 


(c) 


Fig. 2.5. Exemplo 2.2 


A equação de movimento & obtida como se segue: 


aê B A 
EMg=D ”. myL+c(T v)B+k(p v)A+Pv = FD (2.14) 


o mise (BJ Wa Ck(Ã) + Pv = Pr(t) (2.15) 


Nota-se que o fator P/L, na equação (2.15), representa um 
acréscimo na rigidez do sistema e, em relação ao exemplo ante- 
rior, a mudança € devida à contribuição do peso. Neste caso, o 
peso P não pode deixar de ser considerado na equação de movimen 


tn; 


Exemplo 2.3. Invertendo-se a posiçao da viga, fig. 2.6, a equa- 


ção do movimento passa a ser a seguinte: 


ny + c(B) v + (kB) = Bjv = Dect) (2.16) 


Observa-se que a rigidez diminui, em relação ao exemplo 


Ze ls de P/ls 


R ro 
LA 
pare 


Para o caso “em que V=0 e v=0 a equação (2.16) torna-se: 


Pp 


Cedo = Bv = Dr(t) (2.17) 


rio 


Três casos são analisados abaixo, com referência ao pro- 


blema de instabilidade elastica: 


se k > TE , o equilíbrio & estável 
se k < E , o equilibrio é instável, 
se k = Lao , o equilibrio é indiferente. (O 


A carga critica é dada por k AZ/L./ 


O O 
a | o 


| EQUILIBRIO] ESTATICO 


VIGA RÍGIDA 
SEM MASSA 


a 


Fig. 2-6. Exemplo 2.3 
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Exemplo 2.4. Ve-se, pela fig. 2.7, o mesmo problema do exemplo 
2.1, com o acrescimo de uma massa distribuída, m, massa por uni 


dade de comprimento, ao longo da viga. 


VIGA RisiDA 


/ 


(tb) 
p= mg P = mg 
EQUILÍBRIO ESTÁTICO 


tc) 


Fig. 2.7. Exemplo 2.4 


A equação de movimento e obtida utilizando-se o princi - 


pio de D'Alembert, como se segue: 


2Mg=0 7. myL+ U mL VL+k A Av+c L Bv = F(tJD (2.18) 
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(m+ 1 mL)vec(B) vek(A) - 2 r(t) (2.19) 


A unica diferença, na inclusão da massa distribuída, es- 


ta no aumento das forças de inércia do sistema. Se m=mL, pode- 


ça A 


se bem avaliar a diferença entre a consideração de massa concen 


trada e distribuida. 
Exemplo 2.5. No problema anterior, colocando-se a viga na posi- 
ção vertical, fig. 2.8, a equação de movimento torna-se: 


e 2 2 
(m+ DLyvec(p) VÁ Lk (E) + E (P+ 5 pL)Jv = t Flt) (2.20) 


Observa-se, agora, que a inércia e a rigidez aumentam 


mL/3 e a pL, respectivamente, em relação a equação (2.15). 


Fig. 2.8. Exemplo PARE 


Exemplo 2.6. Invertendo-se a posição da viga do exemplo | ante- 
txemplo 2.6. 

rior de forma semelhante à do 2.3 e, consequentemente, o senti- 
do de pe P, tem-se: 


p= 2 a D 
(me Byvrc(B) VALk(Ã) - 7 (P+ 7 LIV = 7 (2.21) 


Ocorre de novo uma diminuição da rigidez. 


- Exemplo 2.7. E apresentado, na fig. 2.9a, um disco com raio rr, 
massa m, espessura h e massa especifica y. São dados ainda, a 
rigidez e o amortecimento 'a rotação'e a torção aplicada no cen- 


tro do disco. A barra que prende o disco e considerada sem 


massa. 


-—— SEM MASSA 


-ky e cr 


Cb) 


Fig. 2.9. Exemplo 2.7 


Para se estabelecer a equação do movimento e aplicado 
> + 
inicialmente o teorema do momento cinetico, ou seja, dHg/dt=To, 
aa SE 


+ 

sendo H., o momento cinetico e LE o momento das forças exter - 
0 

nas, em relação ao centro do disco 0. À parte (b) da figura 


mostra as variáveis que são utilizadas na integração. Assim: 


é i â sn êm pr 
9 ydv = vôn] p? ds=yô | | p?dp p de 
Ho , po vav o 0 do 


td fd 
braço de 
alavanca 


massa do 
disco 


velocidade 
Vinear 
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<> rt 
Ho = 2nY h 6 É 
a 
E e FE o 2 
o +» onde m=nr2h y (2.22) 
-—» 
dH > 
(o) mr? 
e t = o 6 = Las = 7 ama c78- krO 
» mr? 6 c/ô+ k,0 = T 
E aÃ cid 
á 16 ás cr9+ k7º =T ” [2.23] 


Exemplo 2.8. A obtenção da equação de movimento da viga apresen 
tada na fig. 2.10a, composta de duas barras rigidas interliga- 


das por uma rôtula, é vista a seguir. 


A 


pe: E ii a ( DESLOCAMENTO VIRTUAL ) 
Aa Aja 


te) 


o + 


Fig. 2.10. Exemplo 2.8 


A massa da viga, à esquerda da rotula,e desprezível. 0 


deslocamento horizontal e obtido da seguinte forma (fig. 2.10b): 


2v 6V 
o tJê s = 
T = sv. a 6, - E ai (2.24) 
8,= 6 .*. 6h = 26, = 8 (2.25) 


Os calculos das forças, ilustradas na parte (c), são mos 


trados a seguir: 


/ 
Fe = kv/2 
Fa = cv/2 
(2.26) 
Fp = my/4 7 


Fy = (mV-L/2)/2 = my/4 4 
2 

Agora, utilizando-se o princípio dos trabalhos virtuais, ob 
têm-se: 


V V kv 6v 4v dv o 
Z V Z* Ao 6V+ “+ 2 + 2 z+ N RSA e ida 0 


A pin os Ft) É (2.27) 


Observa-se, então, que em todos os exemplos apresentados 
ha um número grande de situações que levam à mesma equação dife 
rencial do oscilador simples, ou seja Ma V+Ca pv +kaçV=F ar sendo Mofo 
Cerf» kof e Foro respectivamente, à massa, Os coeficientes de 
amortecimento e de rigidez e a força externa, efetivos. 
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2.3. SOLUÇÃO DA EQUAÇÃO DE MOVIMENTO PARA SIGL 


2.5.1, VIBRAÇÃO LIVRE SEM AMORTECIMENTO. 


F(t)- o RE 


OU 


2.3.1.1, EQUAÇÃO DE MOVIMENTO 


A equação diferencial a ser estudada, considerando-se nes 


te caso que o coeficiente de amortecimento e a força externa são 
ae SM CENTRAIS os 


nulos, e determinada por condições iniciais de deslocamento e ve 


Tocidade e escreve-se: 


my + kv = 0 (2.28) 
Vv + É v=0 
V + vo = 0 3 (2.29) 


LA - encia angular em rad/seg. 
onde wu, TE E uma requénara SMA 


As condições de contorno em t=0 são v(0)=v, e V(0)=vo- 


2.3.1.2, SoLução DA EQUAÇÃO DE MOVIMENTO 


O processo de solução da equação diferencial (2.29) & for 


necido abaixo: 


A solução geral é obtida por uma combinação linear das 
duas funções dadas na expressão (2.30), ou seja: 
im, t -“ioçt Pula 


(o) ) Ta 
V = c,e 0 + C,e Re: é 


(2.31) 


v- cj(cos w tri sen ngt)+ca cos ugt-i sen ut) 
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v=(c;+c, )cos wgt+(c;-c,)i sen w,t, ou de outra maneira, fazendo 
se A=c,+c, e B=i(c,-c,), tem-se: 


] 
lv=A coswt + Bsenuw é? (2.32) 
a + 


| 
Utilizando-se as condições iniciais, apresentadas no item 


anterior para v=0, obtêm-se A=v e B=v,/0, e assim: 


j———————— = Rá 
/ Va ) 
ny = vo COS ut + To sen uçt Ç (2.33) 


2.3.1.3, CARACTERÍSTICAS DO MOVIMENTO 


CE 


Na solução (2.33), os termos em seno e cosseno são harmô 
nicos com mesma frequência “gp em rad/seg, período To=2n/05' e 
frequência, f sem ciclos/ses. Como exemplo, para ideia de ordem de 
grandeza em sistemas de engenharia civil, indica-se para estru- 
turas rígidas, Tç=0.05 seg, correspondendo a f =20 c.p.s , e pa 


ra as flexíveis, T,=2-00 seg, com f,=0.5 c.p.s. 


Vo cosugt 


(o) 


Fig. 2.11. Composição de harmônicos 


Visto que os dois harmônicos têm a mesma frequência, & 


conveniente compô-los em um harmônico resultante, fig. 2.11-(a) 
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e (b), com a fraguêmeia (8 isto e: 
Et NO 


Dr o O 
2 Fis ” 
vt) =Vortão) cos(wçt-a) (2.34) 
a = tg MS) o detasa Do (2.35) 
Voo 


Pode-se observar que as condições iniciais influenciam so 
E TD SPRESSAO 


mente a amplitude do harmônico. A frequência OE dependente de 


O 


m e de k, sendo portanto uma (Característica do sistema, Ela é 


chamada de frequência natural circular do sistema, sendo (Fo) e 


o» a frequência. eo período naturais. respectivamente. 


2.3.1.4. DETERMINAÇÃO DA FREQUÊNCIA NATURAL PARA SIGL 


A frequência natural pode ser obtida, ou por processo de 


calculo, com a inspeção. da equação de movimento, ou experimen - 


talmente, atraves de instrumentos registradores da resposta, de 


frequencinetros e de analisadores de espectro, Consideram-se, a 


seguir, alguns exemplos associados ao primeiro caso, 


Exemplo 2.9. Para sistema com massa simples, fig. 2.12, a formu 


la, wp=/k/m, pode ser escrita como: 


< 1 a / Dm quê (2.36) 
“o * VYmg P/k Se 
e assim: 
d 0.5 / 
Es ; Lift O 2.37) 
E - w5/2m ' ) o De ) 
ns 


onde 4 É o deslocamento estatico provocado pelo peso P. 
' amento 
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Fig. 2.12. Sistema de massa simples 


Utilizando-se as expressões (2.36) e (2.37), são calcula 
das as frequencias naturais, para dados valores de deslocamento 


estatico, como mostrado a seguir. 


[2090 so 
de = 1/2000 m - f; = 20 c.p.s. (muito alta) 
te = 1/100 m « Ego 8 Cabos 
62 =1m - fç= 0.5 c.p.s. (muito baixa) 


Nota-se que quanto mais rigido e o sistema, maior € a sua 
bn > . e 


frequência. 


Exemplo 2.10. Para calculo da frequência natural, considerando o 
caso de um sistema qualquer, o exemplo 2.8 e utilizado. Da equa 


ção (2.27), tira-se: 


3 ko TON, 
ag = j——=— e dd é 7 m (1+ “EL? (2:38) 
74 m 


Neste caso, a frequência depende também da força horizon 
tal, NM. 

Quando o valor da parcela 16N/kL tende ao infinito, o 
mesmo fato acontece com Usos sendo o sistema infinitamente rigi- 


do, e, consequentemente, o período fundamental torna-se nulo, 


significando que não hã deslocamento na viga. No caso em que à 


mesma parcela e igual a (-1), o valor de vg e nulo, ocorrendo o 


inverso, significando que o sistema levara um tempo infinito pa 


ra retornar à posição de equilibrio, ou seja, e infinitamente 


flexivel, ocorrendo a flambagem do mesmo. Neste caso, então, po 
dd de 


e Ed 


de-se calcular o valor da carga critica do seguinte modo: 


16N 
ee e + As Nor = - TE (2.39) 
, . 
6 k N 
à sk q] 
a o ar 
. vo = [1-N/N...] , (2.40) 


a 


LIMITE ELÁSTICO. 


2 mo): 


Éig. 2.13 processo experimental de Southwell 


Nota-se que wç 2) uma função linear de 1) O processo expe 
ota- 


rimenta] de Southwell, baseia-Se nessa idéia. E 


de (lanbagem com ensaios não destrutivos No regime elastico do 
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. , . f 
material, para diversos valores da carga axial, NAM = re- 


“o € por 


. 2 
gistrando-se os valores correspondentes de tw? 5 0) ê 
01 02º 


fim, atraves de uma extrapolação linear, o valor da carga criti 


cá, Ros e encontrado, fig. 2.13. 


Outros exemplos de calculo de vo são apresentados a se- 


guir. 


Exemplo 2.11. Para a viga flexivel com massa desprezível, figu- 


ra 2.14a, segue-se: 


3... 3EI 
e-Do Seg e ks 
a, = VER (2.41) 


Exemplo 2.12. Da fig. 2.14b, com molas em paralelo, tem-se: 


E. tur 9, ô, - O 


l 


F “ k =k, + Ke > (2.42) 


É a (2.43) 
Us = — à. é bad Pa mo A E 


para as molas em série, da fig. 2.14c, o calculo 


Exemplo 2,13. n | o. 
a - ). dad, FPF.” e6 


e feito da seguinte maneira: 


k,k 
kytk> s" 
= Za 
ser Fr " keq ki+k,> (2.44) 
1 2 


rigidez equivalente para molas em serte. 


25 


142 
vo = (TR (2.45) 


Exemplo 2.14. Para o sistema da figura 2.14d E obtida, a seguir, 


a rigidez equivalente: 


(2.46) 


(c) 


(d) 


Fig. 2.14. Exemplos para calculo da frequência 6 


2,3,2, VIBRAÇÃO LIVRE COM AMORTECIMENTO 


2,3,2,1, AMORTECIMENTO VISCOSO 


a) Equação do Movimento 


No item 2.3.1 foi analisado o comportamento de um oscila 
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dor simples sem amortecimento, fig. 2.15a, submetido a condi- 
ções iniciais, cuja solução & dada pela expressão (2.34). A pre 


sença do amortecedor muda as caracteristicas do movimento, pas- 


sando-se a ter um “movimento harmonico amortecido", ou atê sem 
carater oscilatorio (fig. 2.15-(b) e (c)). 


(a) 


Fig. 2.15. Movimentos: (a) não-amortecido; (b) "harmônico amor- 
tecido": (c) sem carater oscilatôrio 


A equação a ser estudada é a seguinte: 


(2.47) 


35|=m 
< 
) 
fes] 


de E us 
V+=vw 
m » 


Sendo og=k/m e 2Ew,=C/m, vem: 


El 


.. . 2 
V + Eu + UG = 0 (2.48) 


O movimento, agora, é definido pelos parâmetros do siste 


ma w, e E, sendo esta Ultima uma constant. imensional, cujo 


o 


significado é visto posteriormente. 


b) Solução da Equação de Movimento 


A função yaert e novamente escolhida para a solução da 
equação. Assim, substituindo-a na equação (2.48), vem: 
2 


e"tLres2Eu rruç] O rº+2Euçrtu, = 0 .'. 


rr,2 = W9L-EtVEZ-T] (2.49) 
A solução geral & escrita da seguinte maneira: 


2. 4 té k (2.50) 


onde c, e c, São calculadas a partir das condições iniciais do 


problema. 


c) Características do Movimento 


As características do movimento representado pela solu - 


ção (2.49) dependem de três casos, vistos a seguir. 


19 caso - E<l: É o caso das estruturas usuais, correspondendo a 
dm Guav e Sos 


um baixo valor de amortecimento. O sistema é "subamortecido”". 


Então, é verificado na expressão (2.49), que as raízes 


são imaginárias e assim: 
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Es, w (-EtiVT-E7) = -Ev + iu ' (2.51) 
onde 
EA = nar T=E" ) (2.52) 


Substituindo a expressão (2.51) em (2.50) vem: 


-Ew t im. t -iw.t 
v(t) =e º lee * dee *) ava 
Ao. -Euçt o 4 
v(t) = e (A cos w,t+B sen w t) ( (2.53) 
|] ra 


A solução (2.53) representa um movimento oscilatório com 


a amplitude diminuindo segundo a exponencial de (-Ew t) - tem-se 
ne ind do si 
um "movimento harmonico amortecido". Entretanto pode ser notado, 
movimento harmonico amortecido. 
nto, num sentido m 
pela figura 2.16b, que este moviner >, num sentido mais amplo, 
não é periódico. pelo fato da sua amplitude. não ser, preservada. 
Pode-se ainda denomina-lo de "movimento pseudo- harmônico" ou 
"pseudo-periôdico". 
on e q i “ 
Definida na expressão (2.52), w, e conhecida como "fre - 
quência circular amortecida" ou "pseudo-frequência", sendo a re 
WU 
lação 4 « o 
“o 
As constantes de integração A e B são determinadas pelas 


Condições iniciais do movimento, vo e Vo» em t=0, e dessa forma: 


v=v(0)=v => v, = À 


V= V(0) = v, —> vo =-AEuç+Bu, +". B= 


Então: 
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Vo)sen vt] | (2.54) 


mad) 


Compondo-se os harmônicos da expressão (2.54) (figura 2. 


16), a solução pode ser expressa em sua forma final como: 


+ 


/ -Ew Lt Y 
( (0) 2 0 E e 
) v=e gs + = Vig) cos(w t-a) (2.55) 


onde 


(2.56) 


Para o caso de estruturas em que E<l, são exemplificadas 


as seguintes: 


- pontes metalicas E=0.01 a 0.02 
- pontes de concreto E=0.02 
- edifícios E=0.05 

r aço, soldadas E=0.005 


estruturas "off-shore" | aço, aparafusadas  E=0.008 


| concreto E=0.01 


Volta +ENo / Vi -32 


r 
Wê (Vo! ão vo 2")? 


ER) 


Lig. 2.16. Composição do movimento da expressão (2.54) 


30 


20 Caso - E > 7: O radicando, na expressão (2.49) e positivo e 
ad 
as raizes sao reais e negativas. O movimento em questão, nao e& 
oscilatôório e o sistema & dito "superamortecido". Utilizando- 
: a VC 


se a expressão (2.50), segue-se: 


-Eunt wt -wt 9 
V=e [ce  +c;ê J .". sendo e =cosh G+senh O 


e e = cosh 6-senh 6, vem: o 


-Ew t 
v=e " 'cslcosh(u t)+senh (w,t)]+ 


| 
+c, [cosh(w t)-senh(w t)]( 
Para c,+C,=A e c,-C,=B, segue-se: 


v = esugtra cosh(w,t)+B senh(w t)] (E. 57) 


Aplicando-se as condiçoes iniciais em t=0, que são vo € 


—— - a o j 
é vo)senh(o t)] (2.58) 
/Eº-1 


-Ew t 


V 
0 
ja ! Lv, cosh(w, t+ Ei 


A fig. 2.17 apresenta um gráfico da solução (2.58), para 
tres diferentes condições iniciais. Ve-se que em todos os casos 
nao ha vibração e sim um retorno lento a posição de equilibrio, 


devido aq valor elevado do amortecimento. 
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Vo>0 «Voco = 
E SS Vo=0 «Vo>o É 
DS 20" ias V9>0. Vo=0 / 


Fig. 2.17. Grafico de um movimento "superamortecido”, para tres 
condições iniciais distintas 


Os casos da prática não seguem em geral esta situação, 


devido à grande dissipação de energia no amortecimento. 


30 caso - E = 1: Este caso marca a transição entre um movimento 


subamortecido e superamortecido. Na expressão (2.49), as raizes, 


e a solução toma a forma: 


agora, são reais e iguais 
À 


A 
“2 


"q (2.59) 


Aplicando-se as mesmas condições iniciais, ja utilizadas 


nos casos anteriores, segue-se à solução final: 
-ugt , »t 
v=-e EvottVotvoVo ] (2.60) 


representa o limite para movi - 


0 amortecimento critico 


mento nao-periodico e; consequentemente, o movimento retorna ao 


e o o mem me 


repouso, no menor prazo, sem qualquer oscilação, como mostra a 


fig. 2:18. 
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Fig. 2.18. Vibração livre para SIGL com amortecimento crítico 


As partes moveis de muitos medidores eletricos e instru- 
mentos de medida são projetadas empregando-se o amortecimento cri 


tico, para tirar vantagem dessa propriedade. 


d) Significado do Amortecimento Crítico e Coeficiente de Amorte 


cimento 


O significado do amortecimento crítico & mostrado atra- 


ves da seguinte expressão: 


o DaeRO “AR 
= vom | (2.61) 


Quando o valor de E=1, o amortecimento & crítico e 


É é 
Ed Zum, portanto: 


ã ms c pos o RA > up re O id (2.62) 


A expressão (2.62) representa a fração de amortecimento 


em relação ao valor crítico. 


Exemplo 2.15. Fara a equação diferencial, obtida na eq. (2.11), 


do exemplo 2.1, não considerando a força externa F(t), pode ser 


obtido o valor do amortecimento crítico, como é visto abaixo: 


A € vBy". k AS 
V + m (7) Ve + m (7) v=0 


logo, para E=1, vem 


2w m 
[o 2AL 
Ccr = (8)? = “BZ /km (2.63) 
L 


e) Determinação Experimental do Coeficiente de Amortecimento 


A expressão (2.55), para E<l, que & o caso das estrutu- 


ras em geral, foi registrada com à forma seguinte: 
C cos(w t-a) (2.64) 


onde C e « são constantes. 
Da fig. 2.19,que representa O grafico deste movimento os 


cilatôrio, pode-se escrever o seguinte: 


p o o on 
- Fara ut, sm My" 
- Eu t . 
- Para t=t, o come Mad * e o n+1 c e ainda 
+ 
2m x a 
E a EX seudo-periodo) Eai 
tn+l * Es * Ta ua (p P 
2 is E 
Tr a d 
é t + Ste) , 1 
Euo( nova Euçtn VT=E” 
V 1 =C e =( e e á 
n+ 
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“n V1-E2 


n+1 


7 (2.64) 


É interessante observar o rapido decrescimento da rela- 


ção Matt Na com o crescimento de E. Salienta-se que este valor 


não ultrapassa a unidade, no caso em questão. 


Aplicando-se logaritmo neperiano em ambos os membros da 


expressão (2.64), obtem-se: 


V 
Comil ad — (2.65) 


Fig. 2.19. Gráfico da resposta subamortecida - SIGL 


Verifica-se que 6, denominado "decremento: logaritmo", e 


independente do tempo e, portanto, “Util na determinação do va 


lor do coeficiente de amortecimento, E. 


Para baixos valores de E, não e facil medir o de- 


cremento com precisão e, neste caso, pode-se tomar à relação 


Vol/v ' 
nº n+k como: 


Ya Vn Vn41 Vnak-2 , Yn4k-] 


Ynak Yn41 Vn42 = Yn4k-1 Vn+k 
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para k valores consecutivos. 


Como todas as parcelas da multiplicação da expressão an- 


a " | 
terior, são constantes e iguais a 6, entao: - 1º 
v v inata Í 
n 8 ,k n é A V 
=(eº) 2n = k6.:.4€ = 5 tn sam E 
Vn+k “n+k Lu (yr) 
O id 


Pode-se tambem usar a aproximação VT-E? =] na expres 


são (2.65) e assim: 


v 

an —P = 6 = 27E (2.67) 
Yn4] 

Uma outra aproximação é derivada da anterior, ou seja: 


Yn+] 


Desenvolvendo-se expressão anterior em serie de Taylor e 


desprezando-se os termos do 29 grau em diante, vem: « 


“n 
= 1+2785E 


V 
n+1 


“nYn+] (2.68) 


EV n41 


"2 E é 


A fig. 2.20 mostra à relação entre o valor do coeficien | 


te de amortecimento aproximado, com à utilização das expressões 


(2.6)] é (2.68) e o valor exato obtido na expressão (2.65). 


Nota-se que a aproximação (2.67) & bem melhor do que a 


195) 
o 


(2.68) e que praticamente coincide com a exata, nodendo ser usa 
da nos casos práticos de dinâmica estrutural, sem maior esforço 


de calculo. 


E Exato 
E APROXIMADO | 
|O 2038 O SS50 O-9BE0 09806 0.970] 
1.0000 Jp. ——. gg —— go co é e o gp 6 
0 k (1), =----= 
(2) menor 
- = = Le E SD 
00s CIO OIS O2z0o 0.25 E 
exato 


Fig. 2.20. Fator de correção para a coeficientes de amortecimento 
aproximados pelas expressões (2.67)-curva (1) e (2.68) 


curva (2) 


O conrictentu de amortecimento poda atrda ser determinado 
das seguintes maneiras: 
E - amplitude de vibração forçada AM eiméianas dh VOSHOaA Nana 
- defasagem entre a força periódica e a resposta do siste - 
ma ; 


- ciclos de histerese, 


2.3.3, VIBRAÇÃO FORÇADA SEM AMORTECIMENTO 


2.3.3.1, EQUAÇÃO DE Mov IMENTO-SOLUÇÃO 


À fig. 2.21 mostra Os gráficos de uma excitação qualquer, 
30 Yongo do tempo, e o respectivo deslocamento estático, conside 


rendo-se um sistema de massa simples como O da fig. 2.12, 


o 


( b) Vet) = Voos Ft) 


Fig. 2.21. Excitação qualquer F(t) e seu respectivo deslocamen- 
to estático. para O sistema de massa simples 


A equação geral (2.1) € 90 apresentada de forma mais 


adequada, como se segue: 


E (2.69) 
V + 2EuoN + voy à Ea €(t)k/mk 
+ V Ev (2.70) 
N + 2Euçv + voy = wo agf tt) , 


F 
onde . Os 
Rd 


A solução da equação diferencial (2.70) e dada por uma 
soma da solução homogênea, expressão (2.53), e à particular, a 


qual depende da excitação considerada. Entéoé 


(2.71) 


Define-se fator de amplificação instantênea(FAT), ou fa- 


o 
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tor de carga dinamica (FCD), da seguinte forma: 


1 
FAI = LAR / (2. 72) 


a) Carga Retangular 
2.1) Solução Geral 


A equação, para este tipo de problema (fig. 2.22), é apre 


sentada abaixo. 


| flgl= A t>o 
|Elgy-o Elo (2.73) 


2 


.. . z 
V + 2E wav + UV = UG “za 


Velt) 


(b) 


Fig. 2.22. Carga retangular 


A solução & obtida com auxílio da expressão (2.71) e as- 


Sim: 
'y 
-Eu,t 
v=e O [A cos w t+B senwut]+vog” 


one a e a solução particular. 


As condições iniciais são, para t=0, as seguintes: 


v(0) = v,=0 e v(0) =v,=0 


39 
e entao: 
E Em E 
A = Cao e B = E Vas 
e assim: 
-Ew t Ê 
V=v o [l-e º (cos wa t+ é sen vt (2.74) 
bi /T-E7 
carater º 
permanente carater transiente 


A solução geral do problema possui, sempre, resposta nos 


estados permanente e transiente que são responsáveis, respecti 


ts pm 5 


vamente, pelas soluções particular e homogênea. A fig. 2.23 mos 


tra essa caracteristica, considerando-se uma carga retangular. 


Estado 1 
Permonente 


Estado 
E | Transiente | Ê 


' 


Fig. 2.23. Esboço da solução para carga retangular 
Nota-se, pela expressão (2.74), que à resposta transiente & dj- 
“ada pelo coeficiente de amortecimento, E e pela frequencia natural, wo . 
Fazendo-se E=0, na expressão (2.74), tem-se: 


V = VaqLI-cos ug t) (2.75) 
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FAI = []-cos wçt] (2.76) 


A fig. 2.24 mostra o F.A.I. ao longo do tempo. 


FA. 


Fig. 2.24. Fator de amplificação instantânea (F.A.I.) para car- 


ga retangular, E=0 


Se uma carga e aplicada abruptamente, o sistema de amor- 


tecimento não consegue ter influência marcante na resposta. 
Na consideração de E=0, ha o interesse no estudo das car 
9as subitas ou simplesmente pulsos, que são de curta duração,is 


to &, uma parcela pequena do periodo do sistema. 


2.2) Pulso Retangular 


A fig. 2.25 ilustra um pulso retangular com à duração 


t=t,. Dois intervalos devem ser considerados, do seguinte modo: 


Fig. 2.25. Pulso retangular 
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a) Ostst,: A solução & a mesma registrada na expressão (2.76),e 


escrevendo-a de uma outra forma, tem-se: 


is 
Fars T— = [1-cos T tl (2.77) 
eo 


b) t2t,: A vibração é livre e a equação diferencial agora & da- 


da como: 
“e 2 
V + =0 
W 9 


As condições iniciais são as seguintes: 


v(t,) VeoLI-Cos wç5t;] e 


I 
< 
ct 


It) =V, =w,v o Sen ugt; 


A solução da equação diferencial deve ser obtida da se - 


guinte maneira: 


rn yes 


y 
< bord 


: “t 


1 
V=v coswt' + — sen uv t' » para t' = t-t, 
tL, [o o [o 


1 Ca 
V= vagli-cos ut; Jcos vot'+Voo Sen wgt, Sen ut” = 


vt, Ugo t, u ot 
o 
= Veol2sen? — cos ut +2sen S— cos SS sen wt')= 


vota Wgtr 


E Qv, SER (sen ade cos w,t'+cos + sen w t')= 
Vo Lt, (9 t')- 
n HU - 
EV ao sen Sm sen qd 
w t t, 
; E t= -3) 
E 2a SEN —— sen[wç( A J 


nt, t, l 
Eh “sen T senfw,(t- =)] / (2.78) 
0 


mw | 
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max . 1 
= 2sen E | (2.79) 

FA FA/MÁX : 20 

tisTo = 15 


t/To 


ti/Z To - 11/76 


| Lo t/To 


FA FA/MÁX = 2.0 . 
tuZTo = 1.0 FA | FA/MÁX » 2.0 


(d) 


FA FA/MAX * 1.68 FA a 
u/To = MT FA/ MAX » 2.0 


tiu/Z To * 0.5 


Fig. 2.26, Variação do fator de amplificação com (t/To), paraal 


guns pulsos retangulares 


A fig. 2.26 mostra a variação do fator de amplificação 


Com a relação tT, para alguns pulsos representados por ty/To- 


ÀS expressões (2,77), (2.78), (2.79) foram usadas. 


Pode-se observar, nas partes (a) até (d) da figura, que o 


Primeiro valor de Fy/max ocorre durante O pulso e na parte (e), 


q bos, 


“pulso, fm (1) o valor de F/max ocorre no final do pulso. 


43 


A tabela 2.1 mostra os valores de Fy no final do pulso e 


Fiy/RA apos o pulso, variando-se o valor de L/To- 


| V/m 0.5 


Tabela 2.1. Fatores de amplificação para pulso re- 
tangular 


O "espectro de resposta" pode então ser obtido para pul- 


sos retangulares, fig. 2.27. O "espectro" reune parametros de 
ST 


frequencia e resposta maxima para SIGL para um carregamento da- 


do. 


——— 


i 
FA/MAX. F tH 
MAX APOS MAX. DURANTE: E PEQ 0.5 


| MAX. DURANTE MAX. APOS: a <o.5 


Convêr mencionar que como os deslocamentos, forças elas 
ticas e tencies são todas proporcionais, o fator de amplifica - 


c20 pode se; splicado a qualquer dessas grandezas para se obter 
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a relação do efeito dinâmico para o estático 
did , 


Um exemplo interessante para aplicação do assunto ea ob 


tenção do espectro para o pulso triangular, fig. 2.28a. 
b) Carga Senoidal 


O problema e solucionado, a seguir, na mesma sequencia 


em que foi apresentada a carga retangular, ou seja, consideran- 


do-se inicialmente o sistema com amortecimento e uma particula- 


rização para o caso de pulso em seguida, 


A equação diferencial para a carga externa senoidal, mos 


trada na fig.2.28b, e a seguinte: 


Flg) = amwt 
a . 2 1 , 
V + 2Ew,v + Uia E sen wt Ea 
CG 
Tua Zen + ly = wê y sen wt ) (2.80) 
o o o “eo 
Flt) 
(o) 
Fo 


F(t). E sent 


Fig. 2.28. (a) pulso triangular; (b) carga senoidal 
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A solução & dada da mesma forma que a expressão (2.71) » 


onde vy Segue a expressão (2.53) e v E dada por: 


ee: 
- Ç 
Yo * M sen wut + N cos ut / (2.81) 


/ 4 


Substituindo-se a expressão (2.81) na equação (2.80), vem: 


Y p 


a, 
e 
(-w?M sen wt-w?N cos wt)+2Ew (uh cos wt-wN sen wt)+ 
2 gi 
+w (M sen wt+N co = qr 
o s wt) Vo egsen ut 


& 2 VM. E VR 
(ug w2 JM 2EwguN DM ca 


2Ew qorslw ow? )N =0 


é a 1-8? ( 
IM = Veo TI-EERErEE ) S É *"Veo trobtfeeagege | (2.82) 


onde (BB =D) (2.83) 
“o 
be 
A solução global e a seguinte: o ko" 
GS Bens f do e 
a o, oo ad PE a 
TE Mou DD e dE Rega 


pm gp O O 
a. eo E” ; 
v=e (A cos wt+B sen w t)+ T-yjgErpr [(1-B2 )sen wt-2E8 cos wt] 


v (2.84) 


] 
TRANSIENTE y 
SE he 
] 
] 
) 
l 
] 
) 


Fig. 2.29. Esboço da solução para carga senoidal 
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A figura 2.29 


pn a 


mostra a parte transiente, associada à so 
a std E eso 


jução homogenea e a permanente, dada pela solução particular. 


Se o valor de E na expressão (2.84) é nulo, tem-se: 


J” o eee . A 


b.1) Pulso Senoidal -- 


Seja o pulso com meia sengide, fig. 2.30; 


Flt)=Fosenat 


Fig. 2.30. Pulso senoidal 


(2.85) 


Na expressão (2.85), sendo dadas as condições iniciais 


E. v(0)=v,=0 e v(0)=v,=0, para tst,, tem-se: 


AE 1 


seu ia dai dd 
Ni Capeia [sen wt-p sen w. tJ] : 
A NVeo T-B 0 


onde: 


B=D , mas t, = 7m/u ou w=m/t, 
0 a 
o E an 
m/t, 
| ai 
da tm 
AH | 


(2.86) 


(2.87) 
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Substituindo-se a expressão (2.87) em (2.86), obtêm-se: 


, 
o | mt T 

Far —Tora [sen FÉ - qê- sen Erty / GER 

| 1-(geo) Ê 


Para t>t, a vibração e livre e a solução & a seguinte,com 
| | ——— a -, 
fre deita» 


v(t')=c, cos wpt'+c,sen ut! , com V(t=0)=v, e 
V(t'=0)=V, 
Vv=V 2g cos dot cos(w t'-a) (2.89) 
eo B2-T E Gú E <m , 
onde: 


l+cos w t, 
tga = cena 
o 1 


Desenvolvendo-se a expressão (2.89), tem-se: 


É. t, 
E cos(i To) 


- 
H 


A 1? sen[2n(t/T;-t,/275)] » Ou: 
Uitailis 
ho 
a AE 
Fa = FA sen[2n(0/T,-t,/219)] / (2.90) 
onde 
E. JT. enstrt;lt 8 , 
oe e Mescaadhs Fa (2.91) 
1/4-(ty/To) na 
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- Valores maximos do fator de amplificação - Fy/max: 


a) tó 
Da expressão (2.86), com algumas modificações, obtên-se: 


1 t £6 u oda hr 

Fa = T-gr [sen 2mg E -B sen 2y Tl. [ 2. , 

0 valor maximo para o fator de amplificação e obtido, de 
rivando-se a expressão acima em relação a t/T, e igualando-se o 


resultado-a zero, isto E: 


dF 


A 278 t ta. 
FICTA DO) no E Lá (cos 2mg E -Cos 2x E =0 
cos 27B T = COS 27 Te =? 2VB+ ER as 2mt + 2T.m 
0 0 pr To “ 
o | 
t n ) 
. .—— do n=0,1,2,... (2.92) 
, Ta [0] = o SD 
et: 
) too n < ba 
à a 
-1 


(que & uma solução trivial, o ponto na origem, com Fq=0). 


ijto non e E (2.93) 
Oo qu verão ? 
tt; +] 


Os valores de Fp/max sao dados pela substituição da ex- 
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pressão (2.93) em (2.88), obtendo-se: 


= RS ice 

| ID +) 

' F,/max = [o 
A a E do | 

| = (2.94) 


| 
A partir das expressões (2.94) e (2.91) pode ser construi 


da a tabela 2.2, acompanhada das configurações dadas na fig. 2 


31, servindo de auxilio à melhor compreensão do problema. 


MAX IMO 


1.100 
durante 


5.000 mi ad 


2.500 [0.833/1.083 


2.000 |0.800/1.268 


1.083 
durante 


1.268 
durante 


1.500 
durante 


Te TIL 
durante 


1.763 
durante 


0.125 0.317 A 0.493 
apos 


Tabela 2.2. Fatores de amplificação maxima para pulso senoidal 


Pode-se notar que o maximo valor de Fn/max ocorre para 


; apa 


1 - 
7 <1.000 e para n=1, sua determinação e feita como se segue: 


RE. rece 


a) t < t; 


Neste caso fazendo n=1, na expressão (2.94), tem-se: 
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F = L [sen .*. la [o 21 
A max 1 0 y Es a vo (ge; )sen E ] 
ou ainda: 
= 1 2nB 2 | P/n=4 
“A max 1-8 nd Bet A sen Bat (2.95) 


ta) tiZto=- 5,0 


FA 


(c) tiro - 2,0 


(ft) trZTo- 0,750 
FA 
e idda 1.697 


E pe 


| 0.750 1.500 t/To t/To 


FA A (g) ti/70- 0,125 


Fig. 2.91. Variação do fator de amplificação com (q) para pul- 
o 


sos senoidais 


Derivando-se a expressão (2.95) em relação a p e igualan 


do-se a zero o seu resultado, vem: 


dF 
do dyma ir n 2nf 2n 2 2 
BE - q 8º) la COS aiT "Sen ET j ETF cos la 


+ 28[sen a 2m 


A equação acima e resolvida iterativamente, variando-se 


os valores de 8, e o valor encontrado & o seguinte: 


0.625 t, | 
= 0. ou | = 0.8 
E Mesi / Mb 
Do o 
e E q 


Levando-se esse resultado atê a expressão (2.95) obtem - 
se: 


[FpJ/max = 1.768 (maximo maximorum) 


à ta t 


1 


t 
Derivando-se a expressão (2.91) em relação a (po) e igua 
" 


lando-se a zero o resultado tem-se: 


I a sega 
t 
2 - 0.685 | e FK/max = 1.716 | 

To O — 

O valor da maxima ordenada portanto, ocorre durante o pul 
so (t<t,) para - 0.8 com [F,]/max=1.768. 

a O E qe sd 
' 34 


Para qe = 0.5, ou B=1, tem-se + = ! e o valor de FA na 
9 o ' 


n Ou 
0 


expressão (2.95) torna-se , que é uma indeterminação. 


Eira 


qto 2” 
q. tb 
— 


Com a aplicação do teorema de L'Hôpital, supera-se este ti 


po de problema, como & mostrado a seguir 


Ka - 


Tea TÊ A A 

d DO. E 
as + 

TE [sen 278 LE B sen 2y T+ 
pesa pera mes T 

d (1.62) et 

dB B=1 

o tZT 


E 
Fig. 2.32. Pulso senoidal com - = 0.5 
[o 


t 
A fig. 2.32 ilustra o caso do pulso com — = 0.5.Nota-se 
0 DO td 
que o maximo ocorre no fim da meia onda, isto &, o maximo duran 


te e igual ao maximo apos. 


0 "espectro de resposta" obtido & mostrado na fig. 2.33. 


Durante o pulso o sistema comporta- omo rigido e apôs 
t 
como flexivel (a > 0). 
[o 


t, 
No trecho praticamente linear, com O < T— < 0,45. Q sis= 
Nm ça 


tema & muito flexível, ocorrendo uma excitação impulsiva, como 
E o O cánio a O e 


e visto em seguida. 


é - R fd v 
myskv=F , víbração livre com Vo = dv e cuja solução e: v= 00 sen wot (2.96) 


3 
ot) 
[0a] 
— 
| 
O 
ê 
3 
< 
u 
— 
o 
< 
" 
3Im 
o 
+ 
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' 


10 


2º ps nie | TR 
PR A Pd | » 7 To lo 
É ines Na, E: dá E | Es | | 
: / q = — Na — — ê Es | 
| | | | í 
| 


| 
| 
| 


02 os 1 2 3 s 
1%o «fot «1/2 PB 


g. 2.33. Espectro de resposta para força de meia onda senoidal 
SIGL - E=0 adia pan prensa 


Dn o 


Fi 


Anotações 
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Na solução (2.96), Vo=dv para "v=dv, então: 


| dv 

dv=— sen wu t É . Fdt 

| Wo [o dv Mo sen ugt 
No cad, 
0 valor 


do impulso, Fdt, pode ser calculado pela área sob 


o pulso. À fig. 2.34 mostra um pulso genérico e sua área é dada 


por: 


dt = à Fota sendo a uma constante que depende do pulso. Para 


RN Name 


pulsos retangulares, triangulares e com meia senoide, por exem- 


RA, + 
N EA id 


q: Xp + 


plo « tem, respectivamente, os valores de 1, ; e aê. E 


Fit) 


—— Fdt QE, 


Fig. 2.34. Calculo do impulso para pulsos de baixíssima duração 


Utilizando-se do fato anterior, segue-se: 


é Fats a Estao 
a ma E sen uçt dv = é sen wpt 
Fa mn 
sendo 7 *ta é q, 21fg> entao: 

— = de = 2na(fota) 
eo A 
Fa 

TSE * Zn (2.97) 
(6) 2 


[62] 
[Sa] 


(este € o valor da tangente no trecho linear no início do espec 
tro). 


para o pulso de meia onda senoidal, o valor a & a e en- 
as r 


————es em mm 


tão a tangente & igual a 4, que pode ser visto na fig. 2.33. 


Para uma sucessão 


ilimitada de pulsos com meia senoide., 


o espectro de respostas e representado em função do parametro f. 


sua envoltôria & mostrada na fig. 2,35: 


eee 


e = 00/03, 


Fig. 2.35. Variação de Fy/max com 8, para E=0 


Para o caso de meia onda senoidal foi mostrado 


Fy/max(8=1)=n/2. 


Agora, considerando-se m meias ondas senoidais, como na 


/ 


fig. 2.36, para B=1, tem-se: 


Falzm 5 o» m=1,2,... (2.98) 


É interessante notar que o fator de amplificação,na fre- 


quência de ressonância, so tende a infinito se o mesmo aconte - 
» e Se aa aa ni 


cer com o numero de ciclos. Mas e preciso ter em mente que mui- 
tos ciclos podem ocorrer num curto espaço de tempo. Na pratica, 
a ressonância não ocorre integralmente devido ao fato dos siste 


mas não serem completamente lineares. A importante conclusão e 


E a ci AA 


EN 


+ 


que na frequencia de ressonância, ou próxima a ela, as deflexões 
roxima a el: a 


“ - E Lá 
Fig. 2.36 Amplificações para n meias ondas senoidais 


c) generalização do Espectro de Resposta para Pulsos 


Para um pulso genérico como o da fig. 2,39, para tr (tempo 
de crescimento - “rise time") e td (tempo de queda es "decay time"), 
o espectro de resposta pode ser generalizado por um traçado aproxi- 


mado como o incluído na mesma figura. 


O espectro de resposta da Fig. 2.37a mostra que com o aumen- 


to de tr o sistema tende a se comportar como estático com FA=1, ca- 


a. 


so contrário o valor de FA=2 (carga súbita). 
Para pulsos da forma dá FIG. 23763 (tr=t, e td=0), a ten - 
dência & verificar-se o maximo praticamente após o pulso, fig, 2.38. 


Em tais circunstancias a regra tr/T922.5 para [FA tendendo a 


1, fig. 2.39, não Se aplica, pois a mesma & aplicável aos casos em 
que o máximo ocorre durante o pulso. No caso do pulso da fig. 
2.37b1, o aumento de t, o faz tender para a carga súbita, invalidan 
do tambem a regra anterior, pois o espectro tende para 2.0. Então 
ela sô pode ser aplicada para pulsos razoavelmente simétricos, fig. 


2.37b2, escolhendo-se para tr o menor entre tr e td do pulso. 
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(3) 
tr t 
' 
by 
tro 
' s 
º 
Do) ] 2 3 4 s 
ta) 4 /Tos fote a lia ) — e 
. 05 0810 tuna fot 4, 


Fig. 2.37. 


Fig. 2.58. Resposta de pulso em que td=0 


itz 2.5 


Fis. 2.39, Espectro de resposta gencralizado (tr E o menor 


valor entre tr e td) 


Generalização de pulsos o 
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9,3! VIBRAÇÃO FORÇADA com AMORTECIMENTO (PARA CARGA SENOIDAL) 


f 
te, de tal forma que O amortecimento possa ser sentido. Neste 


—— 


caso então, a parcela transiente 


» dada pela solução homogênea, 


desaparece da expressas (2,71) e assim a solução fica sendo: 


p | 
v 
- eo 
TT-87)2+4E7gr [(l-62)sen ut-2ER cos wt] (2.99) 
v 
“eve - cos(wt-a) |, (2.100) 


onde, pela fig.2.40: 


D = [T-B2J2x4E7p7 (2.101) 
& JE 
tg « * TEB. (2.102) 
— 7 “ 
$=7-a (2.103) 


Fig. 2.40. Composição de movimento para a expressão (2.99) 


Quendo o valor de £ é nulo, tem-se: 


0" E qu. 
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à : 
= = T-pr Sen ut . 


F 2.104) 
A “eo ( 


1 
. Fp/max = =.” 


A figura 2.41 mostra o espectro, para o caso de £=D, Ou 


sejas considerando-se a expressão (2.104). Nota-se que quando 
ns que 1 


g=/Z 0 valor de Fy/max e bem próximo a 1, e não 1.24 como obti- 
do para pulsos senoidais. Isso acontece devido ao fato da parce 


ja transiente ser considerada para pulsos, não sendo levada em 


conta para vibração forçada, 


Fig. 2.41. Espectro de resposta para carga senoidal, para E=0 


Na pratica, E=0 não se verifica, mas e Util como indica- 


ção de uma curva limite. 


Vs 
é : tax = 1-8 
. &º 200 Ç =2 
PO] proo ii 
E-= 02 fm 
E= 0.5 2. a -á 


Rm 1, 


Fig. 2.42. Variação do ângulo de fase com os coeficientes de 


amortecimento e frequência 


A figura 2.4? mostra a variação do angulo de fase 4 com 


OS parametros £ e £, utilizando-se as expressões (2.102) e (2.103). 


um 
(o) - 
- -0" a resposta e 
para $=) é TESPOStO esta em fase e para 4=180º em anti-fase. 
o E 


Com à consideração de Ejo, utilizando-se a solução(2.100), 


O E = 


RP ( 
1 
id a: Radios ) (2.105) 


pm 


ET a CD da 
O valor maximo para F tax é obtido para o mínimo de D, 


ou de uma outra forma mais simples 


a 


» Que retrata a mesma conside 
ração, com O minimo de D?, isto e: 


E [(1-82)2+4E2682] = 0 


", -48(1-82)+8E2B = 0 


dar [DE] = -4(1-82)+862+8E2 > O (condição de minimo) 


< O (condição de maximo) 


FA/máx 


Fig. 2.43. Condições de maximo e mínimo para B=0 
Para 8=0 (fig. 2.43): 


der LD] = -4 + 8E? 


Ag > 0+€ > 0.5 = 0.7 + Din *(Fp/max) 


-MBE<0>+E<V/0.5=0.7>D 
max 


para B=0, Ff] 


para B=/T=E7 (fig. 2.44): 


der [D2) = -BE? + 8 - 162 482 =]-2E2>0 


1-2? > 0 + E < V0.5 + Din > (Fy/max) ax 


t=2Eº < 0 E > 5 + não hã (Fp/max) pois 


min 
imaginario 


FA/máx 


—— 


Lo 


p 


Fig. 2.44. Condições de maximo para g=vT-ZgE? 


” (F/Ma) min 


seria 


O espectro pode entao ser construído variando-se os valo 


res de 6 e E, como ilustrado na fig. 2.45. Nota-se que para va- 


lores muito pequenos de E o valor de B tende a 1. 


Quando B=v/T-ZE? tem-se que: 


Fp/max = >—D——— 
TsZESJEsMESUT=ZET]  2EVT-ES 


(2.106) 


Para amortecimentos muito pequenos, pode-se utilizar (o) 


Seguinte: 


+ Sm 
, 
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(2.107) 


FA 


3.0 


o 


FA/mox . 
(ao) SISTEMAS Rígidos 


(b) SISTEMAS SUJEITOS A GRANDES AMPLITU- 
DES E PERIODICIDADE 


E=o fc) sistEMAS FLEXÍVEIS 


Fig. 2.46. Espectro de resposta para projeto 


A nível de projeto, a fig. 2.46 pode ser considerada. 4 


melhor trabalhar na faixa (c), devido à boa participação do amor 
tecimento com a obtenção de uma resposta mais suave e com peque 


nas oscilações. 
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2.3.5» ISOLAMENTO DE VIBRAÇÃO 


Estudam-se agora, dois tipos de problema: A forma como se 


transmite uma força F(t) atraves do sistema para a base; e a que 


regula os efeitos sobre o sistema de um deslocamento da base 


volt). Designa-se para fins didaticos a primeira situação por 


"prapsmiss iba | idades e a segunda por "deslocamento de base",fig. 
qi 


= 


9.41. 


Fig. 2.47. (a) transmissibilidade; (b) deslocamento de 
base 


Caso 1: Transmissibilidade (TR) 
ds a iss 


A força Fb, fig. 2.47a, pode ser obtida do seguinte mo- 
do: 


+ FL = kv + cv / (2.108) 


Um fato importante a ser salientado & que estã sendo con 


Siderado o carãter permanente da força F(t) e portanto, o deslo 


Camento para carga senoidal,F sen wt,pode ser obtido pela expres- 


são (2.100) e assim: 


e. adÚÕÚÇão 


v=- 2 sen(ut-a) = Moo en(ut-d) (2.109) 
D 


substituindo-se as expressões (2.100) e (2.109) em (2.108), 


tem-Se: 
Ts E Nao 
K 
Eu h 'o cos(wt-a) - Cu o 
| E -0) - D sen(wt-a) 
eu . EB 
F p 
. Fo = Cycos(wt-a-y) ê (2.110) 


onde(y é o angulo de defasagem entre a resposta da estrutura e 


o, Suim é 
“tg Y -2E8A ; (2.111) 
———— 
a aan 
ET 
[ea == D ne A > 
E ii 


sendo D dado pela expressão (2.101). 


A transmissibilidade (TR) e definida da seguinte maneira: 


TR = F,/F = Cj = VT+4E2 B? (2.1T2) 
b y max 


A fig. 2.47 mostra a variação de TR com os parâmetros 


e E. TR 
(a) SISTEMAS RÍGIDOS 
3.0 (b) SISTEMAS SUJEITOS A GRANDES 
AMPLITUDES E PERIODICIDADE 
—— (c) sisteMAS FLEXÍVEIS 
2.0 


Fig. 2.48. Transmissibilidade 


utilizando-se do mesmo procedimento que na seção 2.3.4,a 


5 ar ” ; 
abCissa e a ordenada do ponto de maximo ("maximorum") são obti- 


gas pelas expressões abaixo: 


8 =VZ |) (2.113) 


[é 
TR/max = 482 TERRE ora (2.114) 
2A(8E*-4E2-1)+2(1+8E2) 


onde: A = vT+8E 
»6 


Nota-se que para valores acima de /Z o valor da transmis 


sibilidade cresce com o crescimento de E. A melhor faixa de tra 
balho é a (c) pelo mesmo motivo ja apresentado com referência a 


fig. 2.46. 


Caso 2: Deslocamento de Base 
mm E a eç 


Agora, considerando-se a fig. 2.47b, a equação a ser es- 


tudada é a seguinte: e, 
[la 
Ee + FF; + Fa 0 Cal 
Hg ME Re tee 
Pica Ff ? Ju 
“my + cv + ky = 0 6 é 12.115) 


a 


sendo Y=v-v,, deslocamento do sistema em relação à base. 


E ADE 
Nota-se que as forças de amortecimento e rigidez são re- 


e função 


feridas ao deslocamento relativo enquanto a de inércia 


da aceleração absoluta do sistema. 


Desenvolvendo-se a equação (2.115), vem: 


. 2 . 
v+2EugVtugv = ZE Vptudvi f(t) (2,116) 
+ ses 2— .. 
" V+2EugVÁu,V = VA T(t) (2.117) 
f(t). 


onde Vp=“bo 
Comparando-se as expressões (2.117) e (2.70), tem-se a 


seguinte analogia ds Dada pis o o IE] 


2 . —s 2.118 
| bo “eo | 
era ari 

V=v-v, (2.119) 


No caso de E=0 a equação (2.116) torna-se: 


a, 2 2 
vV+uçv = vsYbof (t) (2.120) 


A analogia, fazendo-se a mesma comparação anterior, & da 


da como: 
v v / 
Dt mm. (2.121) 
“bo “eo | 
a 


As expressões (2.118) e (2.121) são muito úteis para so- 


lucionar problemas com deslocamento de base, utilizando os re- 


sultados ja conhecidos de espectros. Todavia, deve-se ter em 
mente que as condições iniciaji do problema em questão, têm que 


Ser as mesmas utilizadas na obtenção de v, quando usada a expres 


são (2.118). 


Para um deslocamento de base aplicado como um pulso re- 


tanigular, uÚtilizândo-sê à analogia (2.121), chega-se à solução 
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jeo espectro de resposta da fig. 2.27 pode ser utiliza- 


2? 
(2 $=0 E o - U- couw, | 
Voa Vro ' 
0. 
É é Cla 
gts [locos vt) (2.122) 
bo 
para um deslocamento de base permanente senoidal, ou se- 
qo onfo * posiel do uma q dale O O 
ia Lá ata MA: o 
Jor 
vp=YboSº” vt + Vp=Vbgu COS ut "a Vbo”-Ybo!” Sen ii 


Utilizando-se a expressão (2.118), se preservadas as con 
a a 


ções iniciais Voz) e Voz), tem-se: 
E aaa RÃ & e 


di 
A co? pet vn Anta he 
ga 
wav / V-v 
o . Msen ut+Ncos ut .'. Em g2[M“sen wt+N'cos wt].”. 
D Vho “bo 


a gr) sen ut+N'cos ut] .-. 


Vbo 
+. prime A) ent)? cos(ut-a")] (2.123) 
Vó B 
Me pr 
tg a! = No 


/ f " 
Substituindo-se os valores de Me N da expressão (2.82), 
X n “ny 
Me 


Nro 


em (2.123), tem-se: 
I 


en : datocam eat 


IA NACW ente 


E. viras Tr cos(wt-a!') = TR cos(wt-a!) par 


= : Ê . ) 
“bo ) Senoio af 


é poge-se utilizar o mesmo espectro da fig. 2.48. 
e utividar O mesmo espectro CE is. Cr. 


Exemplo 2 15. Seja o veiculo da fig. 2.49, que tem uma velocida 
aaa Ses ONA A e 


“e constante, vi, sobre um acidente do terreno de extensão pe- 
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na, 4 com forma senoidal VhoSen ut. 
e , 
qu 


vb = vbo senut 


Fig. 2.49. Exemplo 2.15 


I e o) 


pelo fato da extensão ser pequena pode-se supor que E=0. 


o deslocamento maximo e a velocidade em que o mesmo ocor 
re podem ser calculados com base no espectro da fig. 2.33 e na 


analogia feita em (2.121). Assim: 


tt=0.8 .. LL ft. 0.8 .º. v. = 1.25fd 
1 Vá 


a 
O maximo deslocamento da suspensão da mola deveria ser 
a a 


calculado com uso da analogia (2.118), por ser um deslocamento 


relativo. Mas, neste caso, nao pode ser utilizada devido as con 


dições iniciais A e v. não serem ho êneas como em v CV... 
ESA ses o all o—q— 
Então, é necessário voltar à solução geral e aplicar novas con- 
Ai pasa, 2.85. 


dições a ela. Tem-se: prato 


— bo 2 6 ) 
ug - À cos w,t+B sen uçt- T.gz Sen ut dO (2.124) 


( 


Em pads “o=Yo"pt0) e ué Vostro) 


Vp= Su e. Vo=vo-uvpns-uv, £0 
b7Vbo? €05 ut . “Vo DV bo o É 


ve v » tem-se A- -B " 
com Vo (o) A=0 e B- T=R” USVbo e entao: 


E -8 & 9? v 

nl = T=87 “oVboSen wçt+ T=e sen ut o. 

' a " ( 

- = Fa = Togr [sen ut- g Sen w t] (2.125) 
bo ' 


A expressão (2.125) fornece o deslocamento relativo ocor 


«sendo durante o pulso. O espectro de resposta e visto na fig. 


» 57, notando-se que a ordenada maxima do espectro ocorre apos 


» pulso. Os "maximos apos" são os mesmos do espectro da fig. 2 
à, pera não ocorram com os mesmos valores de ft). O calculo 
sara à obtenção do valor maximo do espectro & obtido da mesma 
:orma da seção 2.3.3.1, item b.1l. Visto isto, o calculo do des- 


tocamento pretendido é fornecido em seguida. 


dis . A PI6 v a TA 


+ velocidade em que ocorre este maximo & calculada da mes 


zm forma anterior, isto e: 
caié E =V. - d 
t,=0.€85 (valor da abcissa do maximQ) .'. vo = 0.685 
1 


v, = 1.46 fd 


k aceleração maxima que o veículo sofre, v , ocorre apos 


EUR 
ruiso € Escime 


e V/ma, vt Vimas "Vilmar = 1.716 u 
A 
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a 


co 


qRECHOS | ——— —— 
NPR o 
fins ur percodiiuidad Map do 
Fig. 2.50. Espectro de resposta para deformação da mola-des- 
locamento de base dy” Vhs senut 


Anotações 


Nota-se que a aceleração maxima tende a zero para siste- 


ivei >0) e a veloci 
exiveis (w OCidade e E - 
cas VE hi Corresponaente E a mesma cal 
iormente Eee co 
curada anterio » COM Oo maximo deslocamento da suspensão. 
- “ 
ã veiculo =u2Vv) & 5 = - 
q aceleração do (v Wav) € igual à pseudo-aceleração do 
ss sed nd 4 fd À 
jstemã. 


pus 


exemplo 2.16. Considera-se, agora, que o veículo do exemplo an- 
Mp ses 


terior tem pelo caminho uma ondulação com varias meias ondas de 


extensão d, ao inves de uma sã e estima-se que o seu sistema de 
amortecedores possui E=0.05. 


teu 
Eu a.. 


0 deslocamento da suspensão e a que velocidade 


ocorre 
são calculados utilizando-se as expressões (2.113) e (2.114). 


e 


mas +" PR = 10.06 Via 


8=0.998 e. fta= 27 e. fz *. v4=1.996 fd 


O deslocamento maximo da suspensão e a velocidade em que 


ocorrem são calculados utilizando-se a analogia (2.118), que po 


de ser usada pelo fato da ação ser permanente, não importando as 


condições iniciais. Assim: 
ceia ni RR AA RA 


2 


wv vi 2 

po = r cos(wt-a) E qá De = E cos(wt-a) 
sa . B2 

e o 


Onde D & dado pela expressão (2.101). 


O "m:ximo maximorum" ou a ordenada maxima do espectro 
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ser calculado da forma Seguinte: 


Fa =0 ou d(F a) 
J a = 


2 
48º L1-82] +4E26º-B"[8E26.4E(1-62)] 0 


D* 
B=0 
o. 4B'(1-B2+282E2) = 0 5 
E me 
VT-ZE? 
para B=0: r 
Sã [A mor ) 
1 a E 
Para B = EDAVO 
/T-ZE ans + 
2 v 4 
Fax = to qo A (2.126) 
B= es bo max 2Ey 1-E 
VT-2E? 
Para valores muito baixos de E a formula seguinte pode 
O din 
ser usada. 
a ] 2.127 
“bo 2 
max 


Para £-0.05 segue-se, então: 


Ee à costas A eine: 1 TELAS sf 
bo o 2:0.05/1-0.0U5* 
en - 1.0025 .'. va=2Bfd 


/T-7=0.057 


v/max=10.0125v,o 


v1=2.0050fd 


Te 
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A aceleração maxima do veiculo & 
AM- 
TR 


va E Tunçe cout -a-4) 
y = -vtvpo TR/max “e V/max = 10.06 ua 
- 2 = 0.998 i aU 
g TA w -998 5 


“. N/max = 10.01 uva DD q 


No projeto de isolamento de vibrações de sistemas & usual 


exprimir-Se O comportamento dos mesmos, em termos de sua "efi 


mm 


ciência de isolamento", designada por (1- 1-TR). Como a faixa dese 


jada para o isolamento da vibração & para B>v2, e ainda na mes- 


ma faixa sendo o amortecimento indesejável, E e possível usar pa 


ra a transmissibilidade a expressão (2.112), com E=0, ou seja: 


1 . 


TR= 2 +". (N-TR) Erê (2.128) 
2-(1-TR) 


2 = 
B T=L TRT » mas 


2 
8? = 2 , sendo va = &( 


onde d, e o deslocamento estático e g a aceleraçã da gravidade 
- MAO O e 


E entao: 
' 2-(1-TR) . a el Te 
“é OS) ud F = ul2m = 7 TR 
É = SET” Qd] é âo [To 2 (TR) 
ai T.UZ5A * de TTOTRJ] .  f=3.1278 de TRI] 
o (2.129) 


Com a expressão (2.129) pode-se montar a carta de projeto de iso 


calculada como se segue: 


«dá K 


vibração, fig. 2.51, 


de 


janentO 


f.nz 


FREQUÊNCIA, 


0.40 0.50 0.60 


DEFLEXÃO ESTATICA , . 
est, in 
gg 


Fig. 2.51. Carta de projeto de isolamento de vibração 
DS LR 


2.3.6. DETERMINAÇÃO DE w, E E COM BASE NO TESTE DE VIBRAÇÃO FOR 
CADA 


Para se submeter um sistema a uma carga senoidal, ja es- 
tudada, a maneira mais conveniente E atraves da utilização de 
duas massas excêntricas, girando a uma velocidade angular w co 
mo mostra a fig. 2.52. As componentes horizontais Ho equilibram 


se, restando as verticais Vo=F; sen ut, sendo Fo=m w?r. 


Aplicada a excitação, registram-se as respostas para se 
ter os valores de 8, atraves de um espectro, como o da fig. 
2.53. 0 valor de B=w/w =] corresponde à ordenada maxima e o va 

õ to 


lor de o e assim expresso pela frequencia w do carregamento, 


s t5 


2.52. Teste de vibração forçada para obtenção de w e E 


Fig. ” 


Bi 1 Bz B 


— tê 


Fig. 2.53. Espectro obtido pelo teste de vibração forçada 


Ymax 1 


O amortecimento deveria ser obtido a partir de = 7Eº 


eo 


mês o valor de vao não pode ser obtido facilmente (B=0?). Entre- 


tanto, para a ordenada v = ie Rag * largura da faixa, B,-B,,» € 


iproximadamente igual a 26. Este fato e mostrado a seguir. 


> v entao: 
(2 /2 eo pia Es" 


V=s = 
2 “max “O ae é 
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e gr - 2(1-262)8? + (1-882) = 0 


nm 


(1-282) + 2EVT+E? 


e para E pequeno tem-se: 


+ 


(1-2E2) da 


T 
o) 
nm 


g? 1-2E-2E? 


mn 


m 
mê 
1] 
Mm 
MY 


B) = 1+2E-2E2 


nm 


t 
cum 
de 
Mo 

aa 


Quando E e pequeno pode-se incluir a parcela E? que o re 


sultado é praticamente o mesmo. Então: 


1-2E+E2 


2 
Bº = E) B, = 1-E 


(148) Bo = 1+E 


nR 


1+2E+E2 


mn 


B> 


Conclui-se então que B,-B,=2E. 


e, cad 


2.4, OUTROS MÉTODOS DE SOLUÇÃO PARA SIGL 


Existenoutros metodos para analise da resposta de SIGL. 


intre eles estão: integral de Duhamel; tecnicas numericas; pla- 


to do ângulo de fase (metodo grafico); tecnicas analogicasstrans 


formada de Laplace. 


São mostrados nas seções que se seguem os metodos da inte- 


ral de Duhame] e uma técnica numérica apresentada por Newmark. 
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n.1. INTEGRAL DE DUHAMEL 


gste metodo esta associado à idêia de impulso, o qual é 


apresentado utilizando-se a solução mostrada na expressão (2. 
54). 
para à fig. 2.544, sendo t) muito pequeno, tem-se: 


Considerando-se agora um sistema que esta inicialmente em 


repouso (v,=V9=0) e que no tempo t=r estã submetido a um impul- 
so I, como mostra a fig. 2.54b, com condições v =0 e Y =Aav=I/m, 

és e 
então a expressão (2.54) torna-se: 


0 , 0 
v E -Ew t* 
: 2] T 0 
v = Pere aid + ) sen w t*]e ' 
sendo t* = t-r 
E -Ewp(t-r) 
Us jÊ Moz sen[wç(t-r)] 
Ce outra forma: 
V= h(t-r)| ; (2.130) 
onde aaa ” 
sen t-7)] -Ew (t-r) 
h(t-+) il o +» I (2.131) 
mu, 


ss 18 


Flt) F(t) 


Fig. 2.54. Impulso I 


A função h(t-r) e denominada “função resposta a impulso 


unitário". Para E=0, segue-se: 
a e Pe 


/ 


 á sen wg(t-t) ( 
ani «e 6-0) (2, 132) 
io ande E 

Considera-se, agora,uma solicitação continua, como ilustrada 


na fig. 2.55. 


Ft) 


Figo 2:55. Solicitação continua 


Para o intervalo infinitesimal, dt, tem-se: 


dl = F(r)dt = Fa f(r)dr, 


se > ane 
ndo Ta variável que define a excitação. 


A resposta elementar e obtida pela expressão: 


”. 19 


Ds Sac E | 
Vá dv = dl h ( t-T ) ( pon Ta aliar (enaça) 
a Z RA 


E. t=t'! 2 
v = , dv = J, F(T)h(t-r)dr 


qet* 
vet), Ve (OM Ct-r)drek Veof UU b(t-r)dr 


Tt * 
E kJ, f(t)h(t-t)dr (2.133) 


Ainda observando a fig. 2.55, se t<t, então t'=t e se 


Bt» Edi 3 


E 


jida a analogia (2.118) e assim: 
ida à analogia À 


- qa kJ f(r)h(t-r)dr 
0 


o : . 
sendo f(t)=VpltT), então: RMS) = (6) 


V = nf. Vplt)h(t-t)dr (2.134) 


E importante salientar, novamente, que as condições ini- 
a O ———<€<*€<<«€.Ç[0C0€0€———0000010———— 


cizis do problema, em v e v, devem ser homogêneas, caso contra- 
no E dE bd E (A Caso 


rio é necessário acrescentar o termo da vibração livre corres - 


ponde v T 
E Vim É Yuar 


xemplo 2.17. Para o mesmo pulso da fig. 2.25, pode-se escrever: 


t 
V = ) vegh(t-r)dr (£=0) 


exemplo 2.18. Para a solicitação da fig. 2.56 
Exemplo Celt 


u 


Us tu 
mM 
t sen w (t-r) 
A 0 lia mo, COsS[w (t-r)] = 
“o ê 
t Se tst, > F.=]-cos t 
coslwç(t-1)]] PA 2o 


0 
se t2t, > Fa-cos[w (t-t,)]-cos w t 


pesenvolvendo-se a expressão para t2t, obtem-se: 
" : 


w t, 


t 
[o 1 
Fi e 2sen —s— sen wo (t- =) 


» OS seguintes re- 


sultados podem ser verificados: 


i)set<t,: o resultado e o mesmo do exemplo anterior. 


ii) se t,stst,: o resultado e o mesmo do exemplo anterior para 


tétsa 


HA) se Tostatss 


t 
1 
Ve IA Foh(t-t)dr+Ih(t-t,) 


iv) se t,stst 


E 
V3 = ve-| Fo Sen ulr-t4) h(t-r)dr 
ts 


Lt, 
ip, ve), o sen wlx-ts) h(t-r)d 


3 


81 


Fosen[w (t- 451] 


Fig. 2.56 » Exemplo 2.18 


2.4.2, MÉTODO NuméRICO DE Nemmark 


Hã diversos métodos numéricos de integração disponíveis 
para a solução da equação de movimento. Os mais usados na prá 
tica são aqueles que Consideram a variação lincar da aceleração 
da massa ao longo do intervalo de inte ção, a qual pode ser 
vista na fig. 2.57, acompanhada das resperira, curvas da velo- 
cidade (quadrática) e do deslocamento (cúbica). 


(E 


* 


Fi : ' 
o 257 Movimento do sistema durante um intervalo de tempo 


(com base na aceleração lincar) 
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A expressão para à velocidade no fim de cada intervalo & 


DR a (2.135) 


para a obtenção de v, em função de v » USa-se a serie de 


caylor» COMO se segue: 


Y o At? o t3 VE. 
vo = vo + VOt+ Vo, Sae + 4 RO é ad de scóti 
Admitindo-se que V varia linearmente então V=(Vo-V,)/At 


avi tal e dal: 


“Vo =V,) + VyAt + É a ai Ve e (2.136) 
As expressões (2.135) e (2.136) são as fórmulas de New- 


No trabalho original, apresentado na "ASCE - STRUCTURAL 


DIVISION - JULHO DE 1959", as expressoes apresentadas são as se 


PS 
tuIntes: 


Va = Vit (DV, Stryvpat á (2,137) 


va = Vatvabt+(S -B)v,At2+ BVo At”, (2.138) 
sendo Eu LU => para reproduzir (2.135) e (2.136). 


ha forma incremental, a equação de movimento & escrita 


mAV + CAVv + kAv = AF (2.139) 


pe (2.135) e (2.136), têm-se: 


% . e! aa RA o At 
Av = Va * Va e tido + AV Em (2.140) 
e 
. .. 2 
dy d Vêm Ma * Vyht + vo A + AV At (2.141) 


Antes de se resolver a equação de movimento na forma in- 


cremental é conveniente explicita-la em Av, assim de (2.141),vem: 
z 6 6 ia 
AV = qe AV qe Ya o 3V; (2.142) 
Agora, de (2.140) e (2.142): 

Av - 3 - SN (2.143) 


Substituindo-se as expressões (2.142) e (2.143) em Et 
139), vem: 
5 6 . 5 3 . AE e E 
lia AV- XE Vj-3V, J+cLe Av-3V,- car Vi J+k Av=AF (2.144) 
Finalmente, carregando-se os termos associados com condi 


des jã conhecidas para o lado direito da equação (2.144), tem- 
se: 


(2.145) 


Onde. 
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| 6 de. 
Femr o tee +) (2.146) 
e sn sta, 
ag 
6. E > 
| AF = aFam(cE V4+3V, )+c(3V, + o V4) (2.147) 


pe posse do valor de Av em (2.145), com a ajuda das expres- 


cães (2.142) e (2.143) obtêm-se: Va=ViHAVS Vo=V,+AV e 
à cjnthv. Este procedimento deve ser feito atê o Ultimo interva 
2 se 


jo de tempo considerado para a integração. 


É importante ressaltar que dentro de cada intervalo At, 
paguraração € Mingar q alem disso o amortecimento e a rigidez 
jo gistena permageCem corstênces, Os erros cometidos são propor 
cionais ao tamanho do intervalo escolhido. Quanto menor o for 
pis prectáão € bbridas Como ele tende a se acumular em cada par 
te do processo, e conveniente no final da integração forçar a 


satisfação da equação de movimento sob a forma integral: 
[E>-F “Pad , (2.148) 
2 


sao as forças de amortecimento e rigidez, respec 


O algoritmo abaixo resume todo o procedimento a ser se- 


guido. 


l. conhecidos v, e v, em At, determina-se F. e Fp, atraves das 
1 1 


relações de força x velocidade e força x deslocamento, res - 


PeCtivamente; 


h/ “ 2y 
Ay pn mta? da .Q ( “401 


v, é calculada em (2.148); 


“ Kear São obtivns, respectivamente, em (2.146) e (2.147); 


>, 
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v em (2.145); 


4 acha-S* ú 
ateu lncam OM Av, respectivamente, em (2.142) e (2.143); 
5. 
' ú Vá + AV e Va = Va + AV; 
6? 
faz-Se Vive e V;=Vz € retorna-se ao passo 1 atê o último in- 
tervalo dt. 
o intervalo de tempo deve ser suficientemente pequeno pa- 
a ser capaz de representar, convenientemente, a excitação, F(t), 


, rigidez» k(t), o amortecimento c(t), e o período fundamental, 


Is 


Normalmente, O intervalo At & estabelecido a partir de To 
pois à representação da excitação não & problema e o intervalo ne 
cessário & facilmente deduzivel e, quanto a k(t) e c(t), o inter 
valo necessário para Ta e, em geral, suficiente para esses ele - 
nentos e pode-se fazer uma subdivisão de At em angulosidades da 
função. Em geral, At=(T,/10) ja produz resultados com uma. razoa- 


vel precisao e, dependendo do caso, pode-se utilizar ate 


st=(T,/50). 


0 método da variação linear da aceleração so & estavel sob 
estavel sob 


: condição de AtsT;/(mvVT-4B), para 8 = & » OU Seja, Ats0.5513T, - 


co 


normalmente. o tamanho do intervalo necessário à precisão é me- 


nor do que este. 
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À siSTEMA COM DOIS GRAUS DE LIBERDADE (261); SEM AMORTECIMENTO 


|. EQUAÇÃO DE MOVIMENTO 


a, Sistema Com dois graus de Iibavdade E nqueis que ne - 
i duas coordenadas pa : 
assita de 0 HÉCêS para descrever seu movimento. 


seja o sistema da fig. 3.1a, definido pelas coordenadas 


duas massas e sráimã 
ne v;» com Constantes Mm em, e constantes elasti 


cas das molas k) e k,. O sistema estã submetido a duas forças 


aternas que variam com o tempo, Fi(t) e Ee TE)s 


= k, 


IR F;(t) kal vo- 4) 
F(t) + 


= kz ma V2 kalvo-v) 
mz 
E(t) Folt) 


(a) (b) 


Fig. 3.1. (a) S2GL; (b) equilibrio dinamico 


Utilizando-se o princípio de D'Alembert encontra-se a 
equação diferencial de movimento para cada uma das massas. Pois, 
perê cada grau de liberdade ha uma equação diferencial para de- 


finir o movimento. Assim, considerando-se a fig. 3.1b, tem-se: 


MV, + kyv, - ka (vz -V4) - lt) = 0 


m;V, 1 Ke bVaeMs “- Fo (E) ns 0 


cc em e ee 
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pearrumando-se os termos nas equações acima, tem-se: 


nda + Chatkadvs - kava = Fy(t) 
(3.4) 


Voy + kavo = Fo(t) 


Estas são as equações de movimento que definem o sistema 


je 26L- 


para determinar a resposta deste sistema deve-se utili - 


ai ds equações (3.1) de forma simultânea, pois v, e v, apare - 
com nas duas equações e são por isso chamadas acopladas. Em ca- 


«o contrário, OU seja, quando estas equações podem ser resolvi- 


jas separadamente, elás são ditas desacopladas. 


Como exemplo de equações desacopladas,e visto o siste- 
ma da fig. 3.2a que e definido pelo movimento vertical v e rota 
ão 6, fig. 3.2b. O corpo rigido de massa m estã apoiado sobre 
duas molas de constante elastica, k, e atuam externamente a 


força F(t) e aomomento M(t). As coordenadas em tal caso são in- 


dependentes. 


F(+) 


Mtt) 
E ===essboassõs 
u o | 


Fig. 3.2. S2GL desacoplado 


As equações de movimento são: 


Ea 88 | 


| 


| 


premplo 3.1. A viga representada DO Pio. $%a 


my + 2kv = F(t) 


0 + 2kd2o = M(t) (3,2) 
e o momento de inercia 
onde Ty ame: da massa em relação ao eixo de 
rotação 


ConstituTda por 


ec o 
duas barras conectadas por uma “otula em B, estã excitada 
por 


quas forças F,(t) e Fa (t) atuando nos Pontos B ec respecti 
, Ctiva- 


mente. O sistema E definido pelos deslocamentos verticais des 


tes pontos e possui uma força normal N at 


Os calculos das forças que estão 


3.3b, estão relacionados à seguir 


re, (c) ! FE» 


Fig. 3,3, Exemplo 3,1 
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myVa » 
E Es = kyvy” 
mo (Vi+V>) 
cm Fes ava 
Vi m,h? Vi, m,h (3.3) 


BT TE Tr 


(Vo -V4) m,h 


E ão Mag 


1,80 momento de inercia de massa de uma barra. As rea 
ções NOS apoios sao calculadas tomando-se os momentos em rela - 
ção aos pontos B. A Ce 5. 


(3.4) 
m 


m 
Ro = (Vavo) - e vy + vo 


N Mm. Mm, 
E (vi-vz) + “E V, + “e Va 


2» 
O) 


Fazendo-se o equilibrio das forças segundo a fig. 3.3c , 


encontra-se a equação diferencial de movimento para o sistema. 


» ct Ro+ e, +Flt).o 


L (m,+m, )V,+ E moVo+(k - Nov,+ É vo = Fylt) 
(é): O (3.5) 


Cart; 


: Mm, V,+ E MoV>+ Ê vat(k>- Adv, = Fo(t) 
Nota-se do exemplo que o sistema & acoplado. 


Exemplo 3.2. A fig. 3.4a representa dois pendulos acoplados por 


Re; ni - 
“io de uma mola, k, que não estã sob tensão quando as duas has 


90 


posiçao vertical. Os deslocamentos angulares anti- 


pstão na 
«nios SãO considerados positivos, 


tes 


porê 


em ue mm 
nº? 


(o) 


Fig. 3.4. Exemplo 3,2 


Tomando-se os momentos das forças em relação aos pontos 


ge suspensão, parte(b) obtem-se a equação de movimento para pe- 


quenas oscilações: 
e 


|] 
o 


ms20, + (mgt+ka2)o, - ka?g,. 
(3.6) 


I 
o 


n420, - ka20, + (mg2+ka?)0, 


acoplada, ; 


3.2, SOLUÇÃO DA EQUAÇÃO DE MOVIMENTO 


3.2,1, VIBRAÇÃO LIVRE 


cmo = e ga 


3.2.1.1, FREQUÊNCIAS NATURAIS E MoDos NATURAIS DE VIBRAÇÃO 


Um sistema com 2GL & caracterizado por duas frequências 
“turais, A cada uma delas estã associada uma configuração do 
“istema segundo a qual este pode oscilar,mantendo-se constante 
é relação entre os deslocamentos, (v, /v, )das massas. Tais configurações 


Cham 


am- : - 
"-se modos de vibração. 
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considerando-se agora o sistema da fig. 3.1a em vibração 


pode-se escrever a partir das equações (3.1), fazendo-se 


pivrê» 
quj=Fa (t)=Do o seguinte: 
pao 
mVa + (kitka)v, - kav, = 0 
(347) 


por analogia com o SIGL o movimento & harmônico e de mes 
' frequência. Os deslocamentos são expressos por: 
, enero: : 


vy=Vy Cos(wpt-a) Vi=-Viuçcos(w t-a) 
(3.8) 
va=Vocos(wt-a) Vo=-Veuçcos(w t-a) 


onde Yy € Ve são às amplitudes características e a o ângulo de 


Substituindo-se as expressões (3.8) nas equações (3.7), 


obtêm-se: 


(-mugtkitka Va ko Vo -=0 
(3.9) 
kaVyt(-mpwd+k> )Vo =0 


Para que haja solução em (3.9), para quaisquer valores 


“hM eVY,, o determinante A tem que ser nulo. Assim: 


(-mwd+ka+k>) ko, 


-k (-mpud+k,) 


Z 


Para facilitar os cálculos na resolução do determinante 


“UpO e : 
P5E-se ky=k,=k e my=m,=m. Encontra-se entao a seguinte equa- 

ro 

E: 
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(8) = ME dO tm O 
as duas raizes são: 
= BT p= 0.382 E 
vo, E G 48) x = 2.618 é 
que são 28 frequências naturais do sistema. 


E 
" 


618 |/K 
o 0.618 Fm 
o. 1.618 VE 


E 
H 


A menor destas duas frequências é chamada de frequência 


a nd a 


sundamental e corresponde ao primeiro modo de vibração. 
mm e 


Substituindo-se os valores das frequencias nas equações 


(3.9), encontram-se os modos de vibração. 


Assim, para vo =0.382 + tem-se a seguinte relação entre 
à 


v k Vo 
he: y = 0.618 e para vo =2.618 — » tem-se: y— =-1.618. 
a 2 


2 


Iº MODO 2º MODO V, J1] 


O Tie “E 
E] “ep (wu) 


so," 0.616 | sos teto [5 K | 14 
õ, = E] “tosa Pe + 4 


1.62 


PESA Modos de vibração €e frequências naturais - S26GL 
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representação dos modos usa-se uma amplitude arbi 


para à 
considerando-se a amplitude do deslocamento da massa Mm; 


; unidade encontram-se os dois modos representados na 


jguê! odem ser chamados d 
= 3.5 e que P e Rh E der 


19 


pro 3-3. admitindo-se, para o pêndulo acoplado do exemplo 
as seguintes soluções: 


o, vycos(wgt-a) 


(3.10) 
vecos(wçpt-a) 


“ 


Oz 


Fazendo-se a substituição na equação de movimento (3.6) 


encontran-Se as frequências naturais e os modos de vibração, que 


580. 


ata Ea] k a? 
Po VÊ ” vo, 2 Vo ta 
VW V, 

Vo V. 


No primeiro modo os dois pendulos tem a mesma amplitude 
ea mola não é estendida; no segundo eles movem-se em oposição 
e aparecera um nô no ponto médio da mola. Eles estão representa 


tos na fig. 3.6. 


Fig. 3.6. Exemplo Basa 
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4. O sistema da fig. 


3. : 
unid | /a possui duas massas diferen- 
a “«aadas por molas, E a ae a 
interliga P le é definido pelos deslocamentos 


ços 


massas e Sua equaça - : 
É vs das quação de movimento, fig. 3.7b, vem a 


ser 
mVa + 2kv4 - kvo = 0 
amv, - kvy + 2kv, = 0 (3. 11) 


utilizando-se as expressões (3.8) como soluções do siste 


aa é substituindo-as na equação (3.11), encontram-se as seguin - 


frequências e modos: 


es 

2 k 2 

Ni OBA vs. = 2.366 

A VV 

ud CL lda V 

r E 16 | qo =-2.732 bo ae) 

v, 

À 


Os modos estão representados na fig. 3.7c. 


Fig. 3.7. Exemplo 3.4 


»” . do) TT 


1.2 RESPOSTA PARA UM S2GL 
plo 


Apos OS calculos das frequências naturais e os modos de 


pração pode-se obter a resposta de um sistema devido às condi 
vi 

5 iniciais. Esta resposta & obtida pela superposição dos mo- 
ae a 

individuais onde cada modo E tido como um STGL independente. 
dos gpa A Da aÃ 
considera-se como exemplo um S2GL sob vibração livre, cu 


- frequências naturais são do, € Vo, e Os dois modos de vibra 


j l = a. 
ção dados pelas amplitudes 4,, e 21» byz e d,,. OS deslocamen- 


tos serão dados pela solução da equação de movimento de um SIGL, 
b da ; 
da bm | 


(8.18) 


E, vibração livre, sendo: 
5 


vi $21€0S(w5 t-a, )tdracos(u, t-a,). 


Vo = $21€0S(w, t-a, +02 Cos(wç tra, dh; 
Desenvolvendo-se os termos em senos e cosenos, tem-se: 


w=0,0,,4C0S w, t+d,4,2 Cos vo, t+*C19,, Sen vo t+ 


“a 1 
+C;0,2 Sen Vo, * 


(5. 13) 


v,=0,02,C0S w, t+d,4,,C0S w t+c, dp, Sen w, t+ 


0; 0, 04 


*C,0,22Sen Vo, E 


Mie dj, do, Cy e c, são constantes a serem determinadas pela 


“Plicação das condições iniciais. 


Em um tempo t=0, o sistema parte do repouso (Vio=V20=0) e 
“ia massa possui um deslocamento inicial V;) € Vo » respectiva 
“nte. Com estas condições substituídas nas expressões (3.13) , 


Enc ) ga 
Ontram-se cy) e c, nulos e d, e d, dados pelas seguintes ex- 


TT 
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TS 
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es 

pr e 
dezVrordaizV>o 
d; 


la PrrV2o-dayVio / 
————————— à dE 


= brrPzz" Przbd> | ii Prrdza dredos ( (314) 


a 


> iii cc Ra 


fazendo-se 921=932=1, tem-se uma Solução geral para d, e 
das à seguir: 
dp? da 


/ Qa2Vao”Vro | d =P22Vio+Vao 
já “ora da Pala = > (3.15) 


Finalmente, com os valores destas Constantes chega-se à 
resposta do sistema que compreende dois harmônicos simples. A 


predominância de um sobre o outro dependera somente das condi- 


ções iniciais. Colocando-se, para o exemplo dado, a resposta sob 


iorma vetorial, tem-se: 


dia biz 
(3.146) 


2 
22 22 


Diz-se que o movimento sô serã periodi quando a rela - 
bo E O em 
ção entre as duas frequências naturais (w, /ws ) for um número 
———— ——— 2 1 ———— — 
O QU 


inteiro, ou se: 


w 
0, n 


i o a . l 
E A, it n; São numeros inteiros. 


Derivando-se as expressões (3.13) em relação ao tempo, 


y =d 
1-0, e +C. uu) Cos w, t+c,uw cos w t 
0,?1Sen Eh devo, Pre sen Vo, é 1 0,911 0, 2 0,12 fo) 


(3.12) 


du 
19 9%) Sen -d, is t+Cjw f,3C05 1 t+C,;w  d,,C0S W t 
0,*21 u E» deu dzz en Vo, 1 0,41 04 20 q, Paz o, 


LES 
< 


considera-se agora que para t=0,2S massas de um sistema partem 


cidades iniciais Vig € Vao, respectivamente, e sem des- 


elo 
on a . e o. 
é ento inicial. Substituindo-se estas condições iniciais nas 
qocaré pe 
vessõeS (3.17) encontram-se d, e d, iguais a zero e ce c; 
exp miiê . 
pads po 
Vij09227Vz0 ( Vistas ar 
( ml - ) as. E 
C4 Vo, b22"P21 | | vg, 0921-622) | 


penplo 3-5+ o pendulo do exemplo 3.2 & posto em movimento com 


seguintes condições iniciais: 0, (0)= A e 0, (0)= 0. As frequên - 


cias naturais e Os modos de vibração Pera cxteutados no exemplo 


3.3. 
partindo-se diretamente das expressões (3.15), encontram 


se os valores para d, e d,. Substituindo-os na equação (3.16) a 


olução é a seguinte: 


A A 
B4(t) = 7 Cos vg,* + 5 Cos o, é 
(3.19) 
A A 
S.(t) = 5 cos 09,8 - 5 Cos o, * 
Reescrevendo-se estas equaçõess,tem-se: 
Wy +W> 
B,(t) = A costa cos(5—)t 
(3.20) 


W +W> 


Ba(t) = -A senta sen(—5—)t 


( 
“omg (w,-u, ) & muito pequeno, 0; (t) e O,(t) tem comportamento 


tua] 0), +U, 0), 4U)2 
à cos( > )t e sen(- gm) y 


Cec 
variando lentamente, segundo mostra a fig. 


respectivamente, com amplitu- 


Toth 


TT 
» 98 


visto que O Sistema é conservativo, existe uma transferên 


cia de energia entre os dois pendulos, 


“premplo 3.6. Considera-se que o sistema do exemplo 3.4 & posto 


em movimento com cada massa possuindo uma velocidade inicial igual 
ato e to 5. Assim pode-se substituir diretamente OS valores correspon 


dentes nas expressões (3.18) e encontra-se: 


c = 0.789 


8) 

nm 
" 

So 

. 

Mo 

mal 

pr] 


À solução & dada por: 


Vão Vio 
V, = 0.789 —- sen 5 t+0.21] sen w t 
ad 1 0, 0, 
Ê , (3.21) 
v v 
Pa TÃO ae gen q teto DL ami a + 
vo, 0; ba é 


em 
Stlo 3,7, Um sistema qualquer de 2GL entra em movimento com 


* Segy ” ã 
SUintes condições iniciais: v,(0)=v,o : Va (0)=v, he Vo = 0. As 
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písticas do sistema são 


ara 
po = e AT — 

= 1.166 Ym nl uv = 1.945 E É PRE 

04 tu, 2 m A fila 

Via 
/ “o / 
1 À ] | E ! Po 
y 2 = “a mn 

a eh -2.414 | 


As equações de movimento são as seguintes: A ) ), 
PRC pa dd vo, t-0.207v, cos 09, + 
(3,28) 


va(t)=0.5Va 0 COS vg, t+0.Sviçcos ug. t 


| partir destas equações pode-se fazer o grafico de Va tEl vs o 
versus To, sabendo-se que Fats fig. 3.9. O mesmo po- 


se-se fazer em relação a v,(t). 


vt 


Vo 
LO 


To) 
W 


.———— e —— e mt 4 


-LO 


Fig. 3.9. 


Ve-se pelo grafico que o movimento & periódico. 
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VIBRAÇÃO FORÇADA 


pd 
esposta de um S2GL 


Ar sob vibração forçada & dada, levan 


, em conta as frequencias naturais e os modos naturais calcu 
; ES fd 
dor 


: e ois à : : 
sos anteriorment » Pois estes são Características do sistema e 
18 


da excitação. 
endem —— 
indeP 


No gráfico da força contra o tempo, fig. 3.10, a carga 


quante em t=r & F(r). Para um S26L excitado por duas forças 

et) e Fo(t), as cargas que atuam em t=r são Fi(T) e F, (7). No- 
ma-se quê area sob a curva em um intervalo de tempo dr é um im - 
ls que causa um acrescimo na velocidade em 1. Assim: 


F, (x) R Fe (rt) 


dv, = EK dt dv, = E Si dr (3.23) 


Fig. 3.10. Impulso F(r)-dr 


Considerando-se, assim, vibrações livres provocadas por 
“es variações de velocidade, eq. (2.32), e fazendo-se a su- 


“rposição dos modos, as expressões para os deslocamentos elemen 


Eres são: 


dv,=c, sen Ugo (t-r)+c,sen vo, Ct-t) (3.24) 
, à 


dvz=C;6,, Sen wo (t-r)+c; dp, Sen vg, (t-1) 
' 
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ec são constantes a serem determinadas. 
pn , “ 
perivando-se às expressões (3.24) em relação ao tempo t, 


ptén-S8 


cw. (t-t)+c;u, cos w (t- 
O 04 0. 0.“ T) 


vw cos w, (t-T)+c,d,2w, cos wy (t- 
gipz6192 1004 0, 0, 04! T) 


pode-se então igualar as expressões (3.25) às (3,23) no 


obtendo-se: 


tempo t=T 
Fa (rt) Fa (rt) Falr) Fs lt) 
m 22 no. m d21- m 
ps si = pe O BE E, 
vg, (922-P21) vg, (921-Pz2) 


Supõe-se que as cargas variam da mesma forma em relação 
E O io o e e —— eee —m -— 


jo tempo, tem-se então: 
Falt) = Faof(T) 5 Falt) = Foof(x) (3.27) 


Substituindo-se os valores das cargas em c, e c,, os des 
locamentos totais de cada massa, em t, é a soma dos efeitos do 


impulso entre zero e t. Colocando-se a solução em forma veto- 


rial, tem-se: 
ER; 


Xá E 
x £1%4) w. f(r)sen w (t-tr)dt+ 
0 1 


o 


FAT 
FAT, 


t . 
+ sda) ,f(t)sen vg, (t-t)dr (3.28) 


(6) 


“nde E e a o ' . , . . - u 
l €. são chamados de "coeficiente de participação modal 


É ué 
“0 dados por: 
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Fi da F>0 Fio F20 
ii És O de m; (3.29) 
= ad oe dl 1 A , 
Er wo, (9227921) vg. (0201-022) 


Eles dependem da distribuição das forças de excitação e 


as propriedades do sistema. 


cada integral nos termos da expressão (3.28) & chamada 

a "fator de amplificação instantânea - (FAI),! e representa a 

iticação de um (S1GI) de frequência ES rtsnão à mesma ex- 
gitação 


A solução final para o deslocamento do sistema & dada a 


cequir: ” : | . 
segu? - ! lana tu se É cfna ) FA 
v = 6,9, (FAID, + Boda (FAI), Ama use | Pp, (3.30) 


; 
fY 


32,2.1. ESPECTRO DE RESPOSTA 


Dado o sistema da fig. 3.11 constituído de duas massas 
r. em, iguais a m, unidas por molas de constantes elasticas k. 


tmessa m, esta submetida a uma força externa, harmônica, 
o 


FSen ut. 


Fig. 3.11. S2GL sob vibração forçada 


E 
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q equação de movimento para as massas torna-se: 


my Va + 2kv, - kv, = Fiçsen ut 
(3.841) 


a + kv, = 0 
mp V2 gi - 


utilizando-se o processo descrito na seção (3.2.1.1), pa 


E geterminação das frequencias naturais e modos de vibração, 


qcontra-S8E 


: om 300 É ini 
= ã — > [o) = 

“o, " 1.62 

E ata 1.00 

w = é - -—+ (o) es 

0, m ú -0.62 


A resposta da equação de movimento segundo a expressão 
(3.30) e a seguinte: 
Va Fio |1.00 Ea 


| 1.00 
- 0.724 (FAI),+0.276 (FAI) 
v. 1.62 E k |.0.62 i 


[3.02) 


O fator de amplificação instantanea & dado pela amplifi- 


cação de um SIGL sob vibração forçada, sem amortecimento. Assim: 
SE A o CO 


(FAI), s— 1 sen ut (FAI), = E sen ut 
(2) - 
a A “os 


Substituindo-se as expressões anteriores na expressão 


(3 
il ton.sê: 
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E 0.724 1 0.276 1 
10 ———» + ———— y|sen ut (3.33 
bles +” 1.172 |1-(-) “Dito. ii 
= (—.) 
vz 91 “o 


a partir desta expressão pode-se construir o espectro de 
+a para o exemplo dado que estã representado na fig. 3.12. 
sta É : a 9. 


se, à partir do grafico, que quando w & aproximadamente 
, na pias 


nota ; 
zero, as duas massas comportam-se estaticamente, ou se- 
qual a == E ia Shi Enc j 
s seus deslocamentos sao iguais ao deslocamento estático da 
io tda % bs 
à, mca 


| 


primeira sa 


quando w torna-se igual a uma das duas frequências natu- 
ais do sistema, às amplitudes das massas tornam-se infinitamen 


te grandes. Verifica-se, portanto, que existem duas condições 


& ressonância, isto e, quando a frequência da força perturbado 
a asc 
ra coincide com uma das frequencias naturais, passando o siste- 


na vibrar no seu correspondente modo de vibração. 


) 
uv. VeIeá 
Ritos 


Cm 


Fig. 3.12. Espectro de resposta para um S2GL 


À partir da segunda condição de ressonância, a frequen- 


rh 
Do = o calt 


TT ————e 
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i ini : ea 
assa à crescer indefinidamente e as amplitudes diminuem 
se nulas . 


AA AS 
putro fator importante que se nota no grafico & quado 


estocomento da primerra massa se anula 4 a frequência da for 
2 : E Em deisés da E 
açitadora torna-se igual a frequencia da segunda mola 

e ) da 
ê 


qeliçTÃo) + Isto quer dizer, que a força atuante sobre a massa 
ur é e- 
“produz somente vibrações na massa m,. A amplitude referente a 


nassa é igual ao deslocamento estático da primeira mola. No 
estã a E a sm — 
cas0» kyckz € É gts dy 

o comportamento descrito acima refere-se a um sistema cha 


ado de "absorvedor de vibrações", em que m, e k, são seleciona 
jos de acordo com a força de excitação e os valores de m, € k;. 


tem-se como exemplo O pêndulo centrifugo. 


0 espectro de resposta analisado é valido para qualquer 
sistema com dois graus de liberdade com F, diferente de zero. 
O ada mm É a sm 


mem 


se, no entanto, estiverem presentes as duas forças de excitação, 


Fe F,, elas deverão ter a mesma função de tempo, f(t). 


ixemplo 3.8. A fig. 3.13a representa um portico cujas vigas são 
rigidas de massas m e 1/2m; a rigidez das barras são k e 1/2k. 
iste sistema & equivalente ao da fig. 3.13b considerando-se es- 
tis mesmas caracteristicas para as massas e as molas. Sua equa- 
ão de movimento & então: 


mv, + 3kv, - kv, = 2F;(t) iai 


mv, - EM + kv, = 0 
Suas frequências naturais e seus modos de vibração são: 


GF id É a k 
0  2m a 2 m 


a. 
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RR 1 0 
(ed A 
f |2.0 Ti 
zm Vo 


pipa 


Fig. 3.13. Exemplo 3.8 


substituindo-se diretamente os valores encontrados ante- 
normente, nas expressões (3.29) e em seguida na solução (3.30) 


tem- Se: 


F F 
alt) = 2 (FAD, + 37 (FAD, 

(3.35) 
valt) = 4 (FAD, = 3 (FAD: 


Vê-se que a contribuição do primeiro modo & cerca de 66% 
iproximadamente do total, pois a influência do, (FAI) e também 


significante. 


demplo 3.9. No exemplo anterior, considera-se que RIM, ESA 


“uma força F,(t) = 4 F,(t) atuando na direção do deslocamento 


b, — e — 
: Sua equação de movimento E: 


mv, + 3kv, - kv, = 2Fy(t) id 


MV. - kyv, + kv, = Fy(t) 


& Cn". 
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o sistema POSSUI àS mesmas caracteristicas calculadas an 


As soluções são: 


arjornente 
F F 
4 10 
alt) = 3% CFAL)s ++ (FAT), 
(3,97) 
Fio “F 
lt) = 3 (FAD, = E (Fan), 


fazendo-se CFAIDD,=(FAI),=1, encontra-se os seguintes des 


qocamentês 
3 Fio 
lt) =7 
F 
15 10 
e 


tomparando-se esses deslocamentos com a contribuição do primei- 
4 10 — . q . .—- 
ro modo (3 -p=)'s ve-se que o primeiro tem uma contribuição bem 


nior devido ao fato de estar mais proximo do modo natural. 


templo 3.10. Considerando-se que o sistema da fig. 3.7 estã ex 
citado por uma força F,9sen wt atuando na primeira massa, sua 


equação de movimento torna-se: 


nv, + 2kv, - kv, = Fyosen ut 
(3.38) 
2mv, - kv, + 2kv, = 0 


À solução final desta equação, considerando-se os coefi- 


ci ira d ” 
'entes de participação e os fatores de amplificação instantaã - 


Tea, e: 
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Fio 0.333 0.333 
va (t) = nóis e, u ) a: w ) sen wt 
A 0, 
(3.39) 
Fio | 0.455 0.122 
já sen 
vel o ça açao a 
0, 0, 
A curva de resposta da frequência estã representada na 
ago 3618 
/ 
vik Vi k 
Fio Fio 
ao ver | | 
LO "o | E 


« 
e sen 


Fig. 3.14. Resposta forçada do Exemplo 3.10 
Comparando-se esta curva com a curva da fig. 3.12, veri- 


lG-se que estes dois sistemas têm comportamentos parecidos, di 


fer . 
“indo apenas nos valores das amplitudes. 
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AÇÃO DE MOVIMENTO SOB FORMA MATRICIAL 
gt 


a equação de movimento, seja para S2GL como para SNGL, PO 


colocada sob forma matricial. Isto visa, principalmente, 
ser 
e 


jet devido a dificuldade de se manusear um grande número de 
$ 


0ó 
pquaçÕeS . 
o cálcul 


esta formulação facilita sobremodo o uso de computado- 


o das vibrações. 
es 
para que Se entenda melhor os diversos conceitos necessa- 


s para à formulação da equação de movimento, sob notação ma- 
rio 

ifal; WET IZA= SE como exemplo um S2GL. Depois estes conceitos 
tn 


- generalizados para SNGL. 
s 


Toma-se como exemplo o sistema da fig. 3.1, cujo equili - 


srio dinâmico para as duas massas pode ser escrito como: 


Fit 
Rg + E. 4 0t) 
(4.1) 
Fr, - FE, a Fa tt) 
Colocando-se sob notação matricial, tem-se: 
E > > ) 
o BR = Etr) o (ia 
ue E 


f - e — . o - 
de Fy e Fe são vetores que representam as forças de inercia e 


: o . - 
lêstica da mola, respectivamente, que atuam no sistema. 


Pela propria definição destas forças obtêm-se a equação 


dá 
“* mov ; : 
“vimento do sistema, sob forma de matriz 


> 
O ak « EUIL (4.3) 
vnde H 

-/ = matriz de massa do sistcehia 
E dá e 

v 


H 


vetor deslocamento 


am 


- matriz de rigidez 
k) si . 


= 


2) = vetor de carga 


A equação (4.3) E a forma matricial das equações (3.1). 


as matrizes apresentadas UMinem completamente as proprie 


do sistema; as matrizes M e K podem ser obtidas 


= a uu Ci, 


des agraniel 
7 apo a 

equilíbrio dinamico das forças, eq. (4.1), ou como mostra- 
0 

pars 


- 


7 MATRIZ DE RIGIDEZ E MATRIZ DE FLEXIBILIDADE 


A força elastica em uma coordenada i, de uma estrutura 


& definida como: 


qualquer» 
= axo Wa 
E, kjj J| (4.4) 
" FE = [7 (4.5) 


nte k;; & o coeficiente de influência de rigidez e é definido 


como a força em i devido a um deslocamento unitário em j e nulos nos 
tais graus de liberdade. , 
A matriz de rigidez pode então ser encontrada por sua pro 


ria definição. Para o sistema da fig. 3.1, fazendo-se v;=1 e 


%=0, fig. 4.1a, as forças para manter estes deslocamentos são: 


(4.6) 


ho contrário, fazendo-se v,=0 e v,=1, fig. 4.1b, as for - 


Gê 
* tornam-se: 


(4.7) 


pode-Se assim, montar a matriz de rigidez do sistema: 


kytko -k, 
(4.8) 


kyxa 


Ee 
eb | 


(b) 


(a) 


Fig. 4.1. S2GL:(a) v,=1, v,=0; (b) v,=0, v,=1 


Outra forma de se calcular a matriz de rigidez e determi- 


êndo-se a matriz de flexibilidade e invertendo-a. A definição do 


“eficiente de influência de flexibilidade, s;;, & o deslocamen- 


vem i devido a uma carga unitaria em j.e nula nas demais direções. 


Consj ça 
'Siderando-se as forças aplicadas Fi .e Fz, os deslocamentos sao : 


é Sj2F, + SyzF> 
(4.9) 
Vs = S21F, + dai 


matricial, tem-se: 


Tm. e 


” . = 
17E 


RR Ri (4.10) 


matriz de flexibilidade. 


2 ——————————— 


E 


nte > 


sqhst Thu indo-se a expressão (4.5) em (4.10), obtêm-se: 


ag kV 
(23. (4.11) 


Kk=1 
di ERES (4.12) 


então, pela definição da inversa de uma matriz e conside- 


ando-se Que as matrizes S e K são quadradas, tem-se da expres - 


s30 (4.12) ua 


S = pesa (4.13) 


Tanto a matriz de flexibilidade como a de rigidez são si- 
útricas e positivas definidas. Estas considerações são justifi- 
atas pelo "Teorema da Reciprocidade" e pela “Energia de Deforma 


ão", respectivamente. 


41,1, ENERGIA DE DEFORMAÇÃO 


A energia de deformação de uma estrutura e o trabalho dos 


sforços em face das deformações, então: 


a 1 ETs 
V=gzzFv;=zFyv (4.14) 


Substituindo-se a expressão (4.10) em (4.14), esta ultima 


“Wrna-se. 


(4.15) 


pazendo-S€ agora a transposição da expressão (4.14), seu 


jo não se altera e tem-se a seguinte expressão: 
gut? 
[ê 


1 STF 
e É (4.16) 


> > 
cabendo-Se que F=Fr e substituindo-se a expressão (4.5) 


(4.16)» obtêm-se para a energia de deformação a seguinte ex- 
pa 


T > 
K v (4.17) 


Finalmente, como se sabe,a energia de deformação & uma 


A Mp e 


3 RARA CAIA 
sencialmente positiva, por isso: 


yrandeza es 


" + 
EEsSF>0 e WKyso (4.18) 


eessimas matrizes S e K sao ditas positivas definidas>pois so as 
sr elas poderão satisfazer as condições anteriores para quais- 


co > 
wer vetores de F e v. 


41,2, TEOREMA DA RECIPROCIDADE 


Supondo-se que uma estrutura esta submetida a dois con - 
nt . 
Os de cargas com seus respectivos deslocamentos resultantes, 
“nsiderand A - . . . 
0O-se que a carga FT e a primeira a atuar seguida da 


“8 Fry, O trabalho de deformação total e o seguinte: 


ER ET Z 


a" 
E «aa RE 


tina parcela não tem o fator 1/2 porque as cargas É, 
otiam com toda a sua intensidade quando ocorreram os des- 
E ex? vd 
pontos VII” 


se agora a aplicação das cargas E invertida, tem-se: 


Fy Fr dt F 


ma 
< 
rm 
ma 
+ 
Doj— 


(4.20) 


os trabalhos W, e W, são iguais pois a deformação da es- 


ra “independe da sequencia do carregamento. Logo: 
tru pe j 


E" 


dE >T > 
' vip = Fr Vi (4.27) 
Este é o Teorema de Betti, que nos diz: "o trabalho das 
nrças do primeiro sistema em presença das deformações causadas 
selo segundo e igual ao trabalho das ações do segundo sistema, 


ssociadas aos deslocamentos causados pelo primeiro”. 


Escrevendo-se a expressão (4.10) para os dois casos de 


arga e substituindo-se em (4.21), obtêm-se: 


*T 


ET > > 
18 fm * Fam (4.22) 


Como os termos dessa igualdade exprimem trabalho pode-se 
fizer 8 transposição de um dos lados da equação, pois assim o 
“alho da estrutura não se altera. Chega-se então à conclusão 
“que Sagl, ou seja, a matriz de flexibilidade & simétrica, 
AsrTE Esta expressão estã de acordo com o Teorema de Maxwell 
tas deformações recíprocas. Fazendo-se estas mesmas aplicações 


Fera 
a : : ut 
matriz de rigidez prova-se a sua simetria. 


dt do 
p vetor das forças de inércia, para um S26L, & definido 
ooo 
F Ma, Mz V, 
a Rs (4.23) 
Fr, Mz1 Mp» V,> 


N- 
ficiente je chamado de coeficiente de influência de 
ndê 0 coe W " 
0 


o E E À s a 77 e 
e é definido como a Éorc) em i devida a uma aceleração uni 
pass pc : 
. a en(á nula nos demais graus de liberdade. 
am iba 
o procedimento mais simples para definir esta matriz e 


seajizar uma estrutura cuja massa total esta concentrada em 
1 A À 


atos cujos deslocamentos de translação são definidos. 


Toma-se como exemplo a viga da fig. 4.2a. A massa de ca- 
a segmento esta concentrada nos nôs 1 e 2. Assim, a massa to- 
:l em cada ponto e a soma das contribuições de cada segmento. 
:tomo somente os deslocamentos verticais são definidos, a ma- 
rizê diagonal. Pois, pela própria definição, a força de iner- 
cia devida a uma aceleração unitaria em 1 & sua propria massa. 


acordo com a fig. 4.2a tem-se: 


Mm, 0 
5 = onde Mm =m + m 
” : la 1b 
0 m, 
Se alem da translação tem-se a rotação como grau de Ti- 
“riade = ' o º 
Cada ponto, pode-se eliminá-lo por condensação esta- 


ita, 


Ma 
sereia Mac 
e 
Mb m 


Fig. 4.2. Massas concentradas em nós de uma viga 


agora, quando a matriz de massa não & diagonal as mas- 


+ são ditas acopladas ea matriz & simetrica e positiva defi- 
$ ao assa ESSE Sa Ee 


nidã. 
gemplo 4.1. O calculo das matrizes de massa e de rigidez, da 


ia do exemplo 3.1, & feito segundo a fig. 4.3 a eb, ce d, 


respectivamente. 


Os coeficientes de massa são: 


ma = (Mmsm,) À 
Ma = Mo e 
2 = m, a 
M22 = Mm, à 
"Os coeficientes de rigidez faz-se um equilibrio de forças 


has 
mola E E 
S Considerando-se a Fig. S.A 


” N7 


eso 
Júj=o 
[em 
N 
=0 
gás 
a | 
h 
(c) o 
N 
feio mem Jvi= é tu + 
| hã 1 fa 
(dg) E 
x nad 
o e 
N 


Fig. 4.3. Exemplo 4.1 


| é | 
| uu |xiz 
| n (1º coso) 
[x | 
h h 
je 2: (2º caso) |xzz à 
7. 


Fig. 4.4. Equilibrio de forças 
= kk. - ki2 = hi 
Kay =" ER h 
N 
kex -r k>22 = k, - h 


As matrizes tem a seguinte forma: 


2(m,+m,) M> 
1 
E mim 
É.» Mo 2m, 
ko. 2N N 
y 1 + h 
á N N 
h ka= | 


Observa-se que as matrizes, K c 4, são simétricas e qua- 
“das 


Tp) : FS E ps 
a, A matriz de rigidez do portico do exemplo 3.8 uv cal 


| 


k 
- ; 3 -1 
a 
k 
z -1 1 
——— LOTE + mio Ko" + 
E ka SE sá Kat 


Fig. 4.5. Exemplo 4.2 


Invertendo-se a matriz de rigidez encontra-se a matriz de 


sJexibilidade. Assim: 


4.3, AUTOVALORES E AUTOVETORES 


Pera um sistema sob vibração livre não amortecida de 
6 a equação de movimento (3.3) na forma matricial torna-se: 


M y + K y = 


1º 


(4.24) 


Colocando-se as expressões (3.4) em forma de matriz, 


120 


(wçt-2) V=-0 wo cos(w t-a) (4, E5) 


cos 


4 o vetor dos deslocamentos modais. Substituindo-se a ex 


- 


n$º a 25) na equação (4.24) encontra-se: 
o (4. 
ess 


, 
Mo +K$=0 (4.26) 


plicando-se os termos dessa equação por M-! e colocando-se 


utti é 
cia obtem-se: 


p em 


(NZK-0GD 4 = 0 (4.27) 


Notando-se que q não pode ser nulo, esta solução so se 


mrna possível se o seu determinante for nulo. Então, 


[MK -oç1] =0 (4.28) 


Esta equação e conhecida como "equação caracteristica" e 


ID 


"suas raizes são os valores característicos, ou autovalores, 


- R N 2 
correspondem ao quadrado das frequencias naturais, no caso, w, 


1 
£ 2 - - 
Eu, . À cada uma destas raizes corresponde um vetor caracteris 
E > E . 


tico é, ou autovetor, que representa o modo de vibração. 
Sabe-se que as matrizes M e K sao simetricas e positivas 
“finidas, mas O produto M-1K=H não apresenta estas caracteris- 


ti Es e : 
is. Contudo, se os autovalores de H são positivos, diz-se que 


Est a EE Ene ' ; EO an — 5 
ê matriz e simetrica e positiva definida; se sao negativos, 
ta a = Siad 
* negativa definida; se alguns são positivos e outros nega- 


tivo — 
“ela E dita indefinida. Finalmente, se todos os autovalo- 


sã e = es a dE 
O positivos e pelo menos um e zero, à matriz e positiva 


“Emj deg: 
“Cefinj = ; — 
Inida; e se todos são negativos e pelo menos um € zero, 


Dm 
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« negativa semi-definida. 


tr? 
nijo; 6 frequencias naturais são sempre valores reais e 


Ê 


C 
é matriz (h é simetrica e positiva definida 


Edo 
post! — 
4.3. considerando-se no exemplo 3.1 que m,=m,=m, k,=k,=k 
apto — o 
o matrizes de massa e rigidez do exemplo 4.1 tornam-se: 
Ne» 
Ê 
1 4m m k 0 
K = 
M 26 m 2m e 0 k 
inversa da matriz de massa e: 
“À ” 9 dl 
4 1 || [O id 
| i im Eh 4h | 


encontra-se a equação ca 


utilizando-se a equação (4.28) 
rcteristica que nos da os autovalores e autovetores, assim: 
of - SE ot + e dr = 0 (4.29) 
vo, = 1.166 VE 6 =: 
0.414 
vg, = 1.945 É de o 
-2.415 
sms 


esta representada na fig. 


À configuração dos modos 


mm 


E Aa, Tp 


Fig. 4.6. Configuração dos modos 
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sTEMAS coM VÁRIOS GRAUS DE LIBERDADE (SVGL) 
$ 


EQUAÇÃO DE MOVIMENTO 


Nos capítulos anteriores tem-se trabalhado apenas com 


e é tido como uma introdução ao SVGL. Neste, as equa- 


o manipuladas de modo semelhante e as dificuldades aumen 
ço 


ã 
tam 

cam: ó número de graus de liberdade. Por isto & que normalmen 
do 


medida que O numero de equações aumenta, ou seja, de acor 


, se trabalha com as equações sob forma matricial. 


por exemplo, para uma estrutura qualquer cujo movimento 
- definido por N deslocamentos, que por equilibrio dinamico en- 
ê 
tre as forças externas e internas resulta em um sistema de JN 


equações, sendo N um número finito. A equação de movimento colo 


cada matricialmente e dada por: 


mM V4 Ni) 
nm A Va) = |Falt) (3) 
ou F(t) (5.2) 


As matrizes M e K foram descritas no capítulo anterior , 


de ordem (NxN), determinadas do mesmo modo que o descrito e apre 
“entam as mesmas propriedades. 
Geralmente, a equação (5.2) é dita acoplada quando a ma- 
tr - : a 
!2 de massa ou rigidez não é diagonal. A diagonalização pode 
Ser = a - 
feita atraves de uma transformação de coordenadas, que e tra 


da pe E 
"O item referente à vibração forçada sem amortecimento, 


Mm 
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Ão LIVRE SEM AMORTECIMENTO 


,, VIBRA 
pur dk uca 
ce um sistema de NGL estã vibrando livremente, pode-se 
É seguintes substituições na equação (5.2): 
; 
qa2é 
ja o,cos (up t-a) o. vm “go, Stem, den) 
F(t) = 0 
2 
Edo, E (5.3) 


Esta equação e equivalente a eq. (4.27) no modo i-êsimo, 


ndo | varia de 1 a N. Ela & verdadeira para qualquer bd; se 
0 
ocorrer O seguinte: 

|n->K-vo, 1] = 0 (5.4) 


Como ja foi dito, esta €& a equação característica onde 
- são obtidos os autovalores e autovetores de M-?kK. Estes, extrain 


a raiz quadrada,são as N frequencias naturais Wo. em Ordem crescente e os + 
A e | e 


tmodos de vibração 4; . 


Assim, o modo de vibração i e uma configuração 4; segun- 


do à qual o sistema pode permanecer em vibração livre, com uma 


frequência w . 
04 


5.2.1, ORTOGONALIDADE 


Supõe-se que uv, eu são soluções difcrentes da expres 


são 1 a 
(5.3) e de 45» seus modos de vibração correspondentes. As 


im 
» Pode-se escrever: 


Mm 
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0; & K 4 (5.5) 


(5.6) 


do-se as expressões (5.5) e (5.6) pela transposta 


escreve-se: 


a T K 
M d; a 05 K 4; (5.7) 
Nos 0,K0 (5.8) 
;M Os 1 É 93 ; 
azendo-Se 3 transposta da expressão (5.7),encontra-se: 
T cs 7 
o, 04; M 1 4d; K 95 (5.9) 
ievido à simetria de K e M. 
Subtraindo-se a expressão (5.9) de (5.8) obtem-se: 
2 2 T : 
oro Edy (5.10) 
mo ug É wo conclui-se que 
1 j 
Ao = é 


“fvidente, então, da expressão (5.7) ou (5.8), que em consequên- 


ci EB 
é da condição (5.11), o seguinte: 


T 
9; 


K +,» 0 (5.12) 


Ex - 
Pressões (5.11) e (5.12) definem a ortoganalidade dos mo- 


ii 
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arã tifo 
5 - a 
P para i-j, à expressão (5.10) sera satisfeita para quais- 
jores finitos das expressões (5.11) e (5.12), que são es 
vã 
” 
qué j 
pitês Guam 
( 
! Md; & Mm; 
(5.13) 
KO 
it É k; são denominados massa e rigidez generalizadas, res 
EPA ; : 


vivamente» associadas ao modo i. 
pé 


outra forma alternativa de se mostrar a ortogonalidade 


mos modos & atraves dos princípios basicos estruturais. 


Admite-se, por exemplo, que um sistema esta vibrando em 


um par de modos r e s, correspondentes a diferentes valores ca- 


- . 2 2 ..= “ 
meterísticos wu, e us. Nas duas condiçoes, considera-se que o 
r s 


úximo deslocamento é atendido e que as massas estão em equili- 


trio estático pela aplicação das forças de inercia, que são: 


E O 1 

(5.14) 
O. M 
e » Ps 


Pelo Teorema de Betti, o trabalho realizado por um deslo 
(am . — . 
"ento virtual correspondente ao modo r quando s esta vibrando 
o 
S4al ao trabalho realizado por um deslocamento no modo s quan 


0 r - ' 
esta vibrando. Assim: 


2 JT 
“o. 0 M Po = vo o” Md 


Ss s 


17 


04H60, =0 (5.15) 
4 uê » quando réfs, então: 


0 (5.16) 


p ortogonalidade em relação à matriz de rigidez & prova- 


ndo O princípio dos trabalhos virtuais. 
u 
” seg 


para um sistema que estã vibrando em um modo s, as for- 


=” inércia e elastica estão em equilíbrio. A força de iner- 
para esta configuração esta definida na segunda das expres- 


TEA a força elastica é a seguinte: 
s0e É 


TA (5.17) 


pelo Princípio dos Trabalhos Virtuais, o trabalho reali- 
ado por tais forças devido a um deslocamento virtual no modo r 


enulo. Então: 


2 A] j 
vg 0 MO, + O K dp = 0 (5.18) 


Ed (5.19) 


42,2, NORMALIZAÇÃO 


Como se sabe, os modos de vibração são obtidos, de forma 


pi RA 
rar q 
là, da equação de frequência (5.3) e por isto não se 


T valor j 
eS absolutos para os elementos de cada velor e sim va- 


Bm RE Es = 
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Pad 


qt 


ava à normalização dos modos pode-se utilizar uma das 


p 


Di ii cada um dos modos por uma constante de 
0 à am 
vo um de seus elementos seja unitario, pode-se utilizar 


pjro termo, mas frequentemente utiliza-se o valor relati- 
im 
pr 


- óxima coordenada; 
ê 
vô 


») este e o processo comumente usado em computadores para 


«iso de vibrações: chama-se &, o vetor normalizado ou modo 
anã 


sa1, obtido segundo a equação: 
porno > 


6; (5.20) 


Tirando-se da expressão (5.20) o valor de 9; e substituin 


na primeira das expressões em (5.13), obtêm-se: 
MB; =1 (5.21) 


use a condição segundo a qual diz-se que os modos estão norma 


lizados, segundo a matriz de massa - ortonormais. 


Reunindo-se todos os modos de vibração em uma matriz qua 
irida em que cada modo constitui uma das colunas, obtêm-se a cha 


Ria matriz modal 6. 


Pia Piz PN 
dz, brz ce. SN 
; à Ló, [o [o ... tu] = Ya Gs; E PN (5.22) 


ET 3 ITA ... “NH 
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T 


cormando-S€  PROSHRES 4 H& u o! Kk 4, os resultados 
qouê diagonalizadas visto que os termos fora da diago- 
tr 
a - go 
o nulOS» devidos às condições de ortogonalidade. Os ter- 


ê io 

1º : E cnaas 

p giagonal correspondem à massa generalizada mp ea rigi- 
â 


d 


Aos yizada ko respectivamente. 
gonerá 


pel ” a 
q reunião dos modos normalizados em uma matriz quadrada 


qa-se matriz modal ponderada. A diagonalização da matriz 
; chá = Ea sa 

| assa por esta matriz resulta na matriz identidade seguinte: 
pe? 


(5.23) 


e, 
= 
ves 
" 
Li 


para a matriz de rigidez toma-se o valor de 4; na expres 


do (5.20) e substitui-se em (5.5) encontrando-se: 


vo M d: e K 4; (5.24) 


remultiplicando-se pela transposta de 6; e aplicando-se a con 


tição (5.21), tem-se: 


ko =u (5.25) 


Utilizando-se a matriz modal ponderada, a matriz de ri- 


ide ' ' 
? torna-se uma matriz diagonal de auto-valores chamada ma- 


“R espectral, isto E: 


* kE q uv? (5.26) 
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o portico de vigas rígidas representado na fig. 5.la 


1. 
sv > as graus de liberdade que definem seu movimento segun- 
st — aunaÃns Vi» Vz2 € Vs. Suas matrizes de massa e rigidez, 
qm, São da forma: 
13º 

m 0 0 2k -k 0 
n=/0 m 0 5 K = -k 2k «k 
Oo O m 0 -k k 


sé frequências naturais e os modos de vibração são deter 


jos pelas equações. (5.3) e (5.4) resultando os seguintes va- 
nã 


jores: 


og, * 0445 Vi vo, = 1.247 JX vo, = 1.802 VE 
1 1 1 

é =| 1.802 de = 0.445 ds = | -1.247 
2.247 -0.802 0.555 


Os Bodos estão ds ode qi na aa 


Ú 
2 
KoF -k Kopt 2k 
t 
2 
k — V| — k,p 2k 0 
2 
Vvo=5 Wood 
2,247 -0, 802 
0,445 


(b) 


0,555 


5.1. Exemplo 5.1 


yum sistema e composto por três corpos rígidos de 


gil» 
18” interVigados por duas molas como na fig. 5.2a. Os cor- 
ms 


ara se moverem segundo as coordenadas apresenta 


| 
pS jivres p 
ão 


as “ira, Sua equação de movimento & escrita como: 
s gigur 
n 


pó 


= 
| - 
ra 
fe) 

< 

a 

o 


0 
0 
n ú i & lek 2Zk <k Vo = 0 (5.27) 
0 
= 0 -k ki Iva 0 


utilizam-se àS equações (5.3) e (5.4) e encontra-se a se- 


te equação caracteristica: 
Mn 
k E 
o = 4 06 * Sir dg * D (8.28) 


0s seguintes autovalorese respectivos autovetores são encon- 


:rados: 


1 1 
vo, = 0 44 — [o = r 
É 
“4 
-1 
2 k 
“ — 3 m [o = 0 
1 


Verifica-se que uma das frequências e nula, tendo para o 


Fodo . - Tas | 
de vibração um movimento de corpo rigido, segundo a 1 TO 
f 


“tb, 
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tb) 


Fig. 5.2. Exemplo 5.2 


«2,3. SOLUÇÃO DA EQUAÇÃO DE MOVIMENTO 


Como ja foi dito anteriormente, o SNGL pode ser definido 
mr N deslocamentos que são associados aos N modos de vibração. 
devido as suas propriedades de ortogonalidade, estes modos des- 
crevem satisfatoriamente, com boa aproximação, os deslocamentos 


com poucos termos. 


De posse das frequências naturais e dos modos de vibra- 
io, a solução geral para um sistema sob vibração livre é a 


“perposição de todos os modos, da forma: 


N 
+ 
Vi E = .º (5.29) 

it) C; bjrostins tem) ê 


un ; ê à 
de à Constante T. representa uma amplitude indelerminada do 
1 


Hovim E 
- 0 1 - , var Es , as 1. 
nto harmônico. Desenvolvendo-se a expressao,tem-se: 
É 0. d.b.cos w,.: 5.30 
Ra Ci;ó;sen a + RA 19; nú 


Ma E e Esse. 
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considerando-se que um sistema possui deslocamento ini- 


+ 
idade inicial v itui - 
» e veloci o” Substituindo-os na expressao 


N 
Vo a ib d;9; 

N (5.31) 
dá Ci un b 
Yo RA 1 0: | 


enltiplicando-se as duas equações pelo produto da transposta 


to modo 9; a qRurIS da MASSA M, então: Bm. 
N 
4 FE 
: (5.32) 
T E A 
Ê i E *d º j=] q Po; 04 


istas equações resultam em uma serie onde os termos que corres- 
cndem a ij são nulos, devido à condição de ortogonalidade dos 
mtos em relação a massa; pode-se encontrar os valores das cons 
“intes que são: 
aee 
- Mov 
e 


(5.33) 
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pa -— 
. termos representam a expansão dos vetores Va ! 
te 


ES 


segun 


do os Seus modos de vibração. 


considera-se que o poôrtico do exemplo 5.1 estã em 


: aci seguintes condições iniciais: 
“” 0 . 


as constantes da solução (5.30) são calculadas pelas ex- 


«es (5. 33), tendo-se c; iguais a zero e para os d, tem-se: 
grosso 


0.543 Vão 


a 
= 
“ 


0.349 vio 


O. 
m 
"W 


d, = 0.108 vão 


A solução final e dada por: 
v(t) = 0.543v, pd, COS wo, t+0.349v, qda COS w, t+ 
+ 0.108v,,óscos w t (5.34) 


03 


demplo 5.4. Para o sistema do exemplo 5.2 considera-se as se- 


vui pé ta 
“intes condições iniciais: 


Fazendo-se as substituições nas expressões das constan 
te 
tes 
“NContra-se: 


Ô 
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d, = Vão d, =d,=0 
Cy = Via C*€, 00 


A solução obtida & dada por: 


+ . 
V = VyodSen v9,t + Vi 09, C0S vg, y) (5.35) 


] 


Vê-se que as amplitudes dependem somente do primeiro mo- 


do. Substituindo-se o modo e a frequência em (5.35), encontra- 


se: 


Vi Vio 
Vo = Vio é (5.36) 
Va Vio 


Conclui-se que o sistema se move, durante todo o tempo, 


com deslocamento e velocidade iniciais. Este é um movimento 


de corpo rigido. 
5.3. VIBRAÇÃO FORÇADA SEM AMORTECIMENTO 
vação de movimento de uma estrutura com 


A solução da eq 
NGL pode ser encontrada de duas maneiras. À primeira segundo a 


Cm 


mesma formulação expandida apresentada para um S2GL; a segunda 


é conhecida por equação modal e e a mais comum nos livros de di 


namica existentes. 


5.3.1, FORMA EXPANDIDA 
Supõe-se que uma estrutura estã submetida a N forças ex- 


ternas, definidas como: 
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26) = Fo Ot) (5.37) 


E parte do repouso. De acordo com a fig. 3.10, encontra-se 

q : : 

' acréscimo de velocidade a seguinte equação: 

rã 

pê 
> E ) 

+ -1 dr/ 

df o MFo fluid (5.38) 

os deslocamentos elementares são dados, considerando-se 

je cada uma das massas esta sob vibração livre depois do acrês 


mo de velocidade, sendo: 
[o] 


N 
+ 
K pa cidjsen E EE) (5.39) 
Derivando-se esta equação em relação ao tempo e igualan- 
do-se à expressão (5.38) em t=t, tem-se: 
N 1Ê 
) Cidjuo. = M Fa f(r)dr (5.40) 
Premultiplicando-se por 63M e aplicando-se as condições de ortogo- 


nlidade para i£j, encontra-se: 


a É 
q 3 “a 

c,= q IO f(r)dr | (5.41) 
“a 95 M 95 | 


Ea . 
Substituindo-se a expressão (5.41) em dv (5.39) e integrando- 


Ç - . . 
“Os membros, encontra-se a solução da equação de movimento. Assim: 


w. f(r)sen w, (t-v)dr 
] Us 4; M d. J, “a ps 


(3.28) nota-se que a uu pa 


| MM 
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o “fator de amplificação instantânea - (FAL) ,º - e é definido do 


mesmo modo que na seção (3.2.2). Então a expressão (5.42) pode 


“ser escrita como: 


> N E / 
V = Cc. 6. / 
sy Ca SUA; ( (5.43) 
onde, o ; 
- 65 É, 
É eg (5.44) 
“504 6% 


e um coeficiente de participação e relaciona trabalho de defor- 
: e a O mm ee, 


mação devido às forças externas com a energia cinetica. 


a 


5.3.2. PROCESSO DA EQUAÇÃO MODAL 


A equação modal é caracterizada por uma transformação de 


coordenadas cuja finalidade é desacoplar a equação de movimento. 


Sabe-se que os modos naturais são vetores linearmente in 
dependentes e que podem definir por uma combinação linear, qual- 
quer vetor do espaço considerado. Assim, o vetor dos deslocamen 


tos pode ser escrito como: 


V = dam, + dano + cce + Oy (5.45) 
+ > 
ou v =0n (5.46) 


onde é & a matriz modal que transforma as coordenadas generali- 
+ = à > 
zadas n nas coordenadas geometricas v. Estas coordenadas genera 


' = (| 
lizadas são chamadas de coordenadas normais. 


Derivando-se a expressão (5.46) duas vezes em relação ao 
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tempo € substituíndo-se na equação de movimento (5.2), encontra 


se: 


Mons+kKone=F(t) (5.47) 


Premultiplicando-se a expressão (5.47) pela transposta de 


um vetor no modo i, torna-se: 
+ 
kom =6; Ft) (5.48) 


Se os termos antes da igualdade forem expandidos, todos- 
os termos diferentes de i desaparecem e o resultado & o seguin- 


te: 


Loo e oLE ama cATÃtO 5.49 


Das expressões (5.13) tem-se: 


ã T+ 


Dividindo-se a expressão (5.50) por m,. encontra-se finalmente 


a equação desacoplada, isto &: 


+ 
94 FC) 
da o Ps Da (5.51) 
ou 
T 
A: + wu? = qu? a Rá. f(t) 
Rai Dado ERP US (5.52) 
Í tugd NO; 


Comparando-se a equação (5.52) com (2.70), para um SIGL 
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e o seguinte: 


qui 
o" 
[or F = tas A , al 
Vo as | 2 pas Cx paceu ) 
à Í (5.53 
2 . 
Neo; Po; 4; M 44 / iai n mo pepandido 


[sto corresponde ao deslocamento estático segundo a coor 


tiram 


en . 
i. 


ação entre Na 


A expressão de Neo, em função da matriz de rigidez pode: 
” conseguida atraves de um artifício, onde se coloca entre o 
ção O produto de KK-*, sem alterar o resultado final. No denomi- 
rador substitui-se o valor correspondente segundo uma das expres 
sões (5.7) ou (5.8), quando i=j, obtendo-se: 


0%; A 
ai GR 9 (5.54) 


“Mo, * | (5.55) 


Conclui-se portanto, que Neo e uma expansão segundo os 
es AR ES 
modos de vibração do deslocamento estatico Vaô - devido ao car- 
"egêmento F, - por meio da matriz de rigidez. 
O deslocamento final & dado pela expressão (5.45), onde 


nê determinado pela equação (5.52). Este pode ser escrito | da 


Seguinte forma: 


(6 RA 
v, = 6; Neo (FA); (5.56) 
Gy —+ N E ' 1) 
V = E Da did 
ty Meo, MFA, 
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quação modal pode-se também trabalhar com a matriz 


Na € 
qual izada» fazendo-se V= õ n. Substituindo na equa - 
8? vimento (5.2) chega-se à equação final: 
«08 
O 
+ Tt 
Feune 2 f(t) (5.58) 


sãa matriz espectral. 
qnde Vo 


para um modo r qualquer, esta equação toma a forma: 


e 
A+ vê àr= Op Fo f(t) (5.59) 


pemplo 5.5. 0 portico do exemplo 5.1 & excitado por um vetor 


e forças» fig. 5.3. A solução da equação de movimento, segundo 


; forma expandida, e feita utilizando-se as expressões (5.43) e 


(5.44). As constantes c; são: 


F 
(8) 
| = 6.164 2 


ol 
" 


a) 
ty 

H 
o 
. 
fem) 
—s 
[6] 


A solução final é a seguinte: 


< 
" 


F F 
6.164 1 44 (FAI),-0.18 E da (FAI),+ 


F 
+ 0.018 “2 q, (FAI), (5.60) 
Com a expressão dos (FAI),. para uma carga senoidal de 
Ig 
(tem-se as três amplitudes. Nota-se por esta equação «que a 


0n : me e S 
tribuição do primeiro modo e maior do que d dos outros dois. 
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postando ss este exemplo pela equação modal, pode-se 


qiretamente à expressão (5.57), onde tem-se para 


» os 
ante intes valores: 
e; F 
- 6.164 E 
Neo, 
a =D: 8. 
Ea 
neo, + 0018 


romparando-Se com as constantes calculadas anteriormente, veri' 


Fica-se que elas tem os mesmos valores. O resultado final sera 


ontão O mesmo que o anterior. 


conclui-se então, que os dois processos são validos 


chega-se ao mesmo resultado. 


——e Vo f(t).sentwot 


Fig: Das Exemplo 5:65 


>. VIBRAÇÃO FORÇADA COM AMORTECIMENTO 


À equação de movimento para um SOL, sob vibração força- 


Con - ds ' 
m amortecimento, é a seguinte: 
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KV = F(t) (5.62) 


- 


ot chamada de matriz de amortecimento e seu 
tm . 
nê 


fluência de amortecimento 6, , sé definido como a 
1 seu 


coeficien 


força 


e : : : en 
; ordenada i devido a uma velocidade unitãria na coordenada 
co 


nê 
ie 
neste sistema, os modos de vibração são os mesmos que pa 
E, 
it sistema não anortecido (dependendo a formação da matriz 
rã e — - 


b. 


cl. CONDIÇÕES DE ORTOGONALIDADE PARA O AMORTECIMENTO 


A condição necessaria para que os modos sejam ortogonais 
- matriz de amortecimento e: 


00800 *0 (5.63) 
sara mén. Para m=n, tem-se: 

Midas (5.64) 
e “n-n n e 


we e o amortecimento generalizado que pela propria 
CT me e 
tn 


definição 
“o coeficiente de amortecimento e: 


; arma tod a S46L 
” ESA E Ju, 
cn = E, é ni vm E (5.65) 


PRayleigh definiu a matriz de amortecimento como: 


/ 
Cc = 
E àQ M + à, K | (5.66) 
e q e a, são fatores de proporcionalidade arbitrarios. 
| 
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Esta equação satisfaz a condição desejada (5.63) apenas mul 
tiplicando-se ambos os lados pela operação de ortogonalidade.Qual- 
“quer matriz que seja uma combinação linear das matrizes de mas- 


sa e rigidez e ortogonal. Assim, a matriz de amortecímento orto 


gonal pode ser escrita da seguinte forma: 


cn] a CM-2KID ; Go) (5.67) 


tm 


O amortecimento de Rayleigh esta naturalmente incluído 
nesta equação. Com este tipo de equação torna-se possível calcu 
lar os coeficientes de influência de amortecimento necessarios 


para que o sistema seja desacoplado com coeficientes de amorte- 


cimentos desejados. 


Para o amortecimento generalizado, em um modo n, a con - 


tribuição do termo b da série, dado pela expressão (5.67) & o 
seguinte: 

Toco =açd! MEM-IKID O (5.69) 
n-n b?n == — n pi 


Da equação (5.5) sabe-se que K é =us M 4. Premultipli- 


n 
cândo-a por o! k M-1, tem-se: 
Ty RD aah 
On K MK On = ug dn É Pn Vo, (5.69) 


Por Operações semelhantes a estas pode-se mostrar o seguinte: 


67 mena 6, = vo? Mm (5.70) 


DE caga 


. 2b N 
. . /Enb = *b don A (5.71) 
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assim, à matriz de amortecimento para algum modo n e dada 


nO? 


2b 
“5 Ca “4 da Mw M 5.72) 
; 1 (2b-1 
ed im E ) “pvo, ) (5.73) 


ps, - 


esta expressão (5.73) proporciona uma maneira de se ava- 
iai 03 CONSTANTES Ap para dar os coeficientes de amortecimento 
jsejado em algum modo especificado; as constantes são determi-. 
mdas então, das equações simultâneas para os diversos coeficien 


ves. À princípio, b varia no intervalo (-e,5), mas na prática 


procura-Se utilizar os valores mais proximos de zero. 
a me E - - ——— 


É interessante notar da expressão (5.73) que quando a ma 
triz de amortecimento e proporcional à matriz de massa (C=a  M; 


»=1), o coeficiente de amortecimento & 


inversamente proporcio - 


nl a frequência de vibração; então os maiores modos de uma es- 
a Pie Aesttiasdes pes a LN ps 


trutura serão pouco amortecidos. Agora, se o amortecimento e 


proporcional a matriz de rigidez (C=a,k; b=1), o coeficiente de 


imortecimento & diretamente proporcional a frequencia; e os 


misres modos da estrutura serão lentamente amortecidos. 


Como exemplo, deseja-se avaliar as constantes para pro - 
Porcionar tres coeficientes de amortecimento específicos. As 


“quações resultantes da expressão (5.73) são: 


é Ro lis Uy E 
UI W1 
E 111 1 (5.74) 
(à = 5 |—=5 — u) a e 
: 2 |w; uz à 0 
Ea e da DE a, 
DE Us Ê: 
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ão DA EQUAÇÃO DE MOVIMENTO 


soLUC 
qu 


perivando-Se a expressão (5.46) uma e duas vezes em rela 


tempo» substituindo-se na equação (5.59) e seguindo-se O 
+. à0 
30 


, 


procedimento do item 5.3.1, chega-se a equação final, is- 


PE 
mm é: 
+ 
e 2 2 o; LH 
je + CE quo. e Po Mi = “o, E Bud f(t) (5.75) 
“o, 44 = d4 
. 2 à 
” hj + Eq, MA + vo Ná = “o Neo, f(t) (5.76) 


Essa equação é equivalente a equação para um SIGL sob vi 
2. EO 


(279). Sua 


alução & dada segundo a expressão (2:82), onde as constantes A 


plo Ss Awt 


«ração com amortecimento para carga senoidal - eq. 


:«: são determinadas pela aplicação das condições iniciais. A 


reSDOST 


tz final serã dada pelo metodo de superposição modal. 


vesplo 5.6. Para a estrutura do exemplo 3.8, define-se uma ma- 


" a 
. 


acortecimento tal que os coeficientes de amortecimento 


são 54 do crítico. 
Dando para b valores entre 0 e 1, e com os valores das 


recuências calculadas naquele exemplo, tem-se pela expressão 


3.74) o seguinte: 


0.05 1.414 VE à, 
| 
a 

0.05 0.707 VE a, 
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a) [0.047 É 


ao jo.087 fl 


A matriz de amortecimento segundo Rayleigh & então: 


0.235 -0.047 
C = 
“2 U-0.047 0.094 
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